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ABSTRACT 


RESEARCH RELATED TO MULTI DEGREE OF FREEDOM 
MAGNETIC SUSPENSIONS 

Yan Yang 

Old Dominion University, 1997 
Director: Dr. Colin P. Britcher 

This thesis covers the activities associated with 1) recommissioning of the 6- inch 
Magnetic Suspension and Balance System (MSBS) and 2) enhancing the Annular 
Suspension and Pointing System (ASPS). These activities continue the program of 
research in the multi degree of freedom magnetic suspension program at Old Do mini on 
University. The 6-inch MSBS is a large gap magnetic suspension system used as a model 
support and balance system in wind tunnels. The first stage of recommissioning of the 6- 
inch MSBS for the 6-inch subsonic wind tunnel is performed. Experimental data and 
computational results for the magnetic field distribution of the MSBS are given and 
compared with the estimated magnetic field strength. The comparison indicates that the 
magnet system is still operating as designed. The Electromagnetic Position Sensor used in 
the 6- inch MSBS is analyzed before its final reinstatement. The ASPS is a small gap 
magnetic suspension system providing orientation, mechanical isolation and fine pointing 
of space payloads. The control system of the ASPS is improved by designing a new code 
to realize real time control over the system. The original Proportional-Derivative 
controller is upgraded to a Proportional-Integral-Derivative controller. This progress leads 
to a system which is more stable and robust with less noise. This new controller is generic 
and can be adapted to the 6-inch MSBS. Finally, future work on these two systems is 
proposed. 
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CHAPTER I 
INTRODUCTION 

LI Magnetic Suspension and Levitation 

The capability to generate non-contact forces has attracted scientists working in 
many different fields. Practical magnetic suspension systems first appeared in the early 
1950’s. Over the last several decades, developments in control theory, innovations in 
sensing instrumentation, and advancements in digital computer hardware and software 
technologies have all helped to bring this promising concept to industrial application. 
Nowadays, magnetic suspensions are used for bearings of gyroscopes, as plasma 
containers in proposed nuclear fusion electric power generators, as supports for high- 
speed ground transportation systems, as interference-free model support systems for use 
with wind tunnels, and many other applications. 

Long-term efforts have been made to utilize magnetic suspension concepts in 
aeronautical and astronautical applications. A series of research projects sponsored by 
NASA has been conducted since the 1960’s [1]. These applications include both small and 
large gap suspension systems. Small gap magnetic suspension system activities at the 
NASA Langley Research Center (LaRC) began with the concept of the Annular 
Momentum Control Device (AMCD) and developed into the Annular Suspension and 
Pointing System (ASPS), which is designed to provide orientation, mechanical isolation 
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and fine pointing of space payloads. Large gap magnetic suspension systems have mainly 
been developed for use in wind tunnel model support systems. The research at LaRC 
began with the development of a single degree of freedom demonstration system. Since 
then, several Magnetic Suspension and Balance Systems (MSBS) such as the 6-inch and 
13-inch models have been developed under various grants and contracts. Future 
requirements for systems of this type are associated with test activities such as elimination 
of aerodynamic support interference in wind tunnel testing, gravity simulation of store 
jettison testing, etc. [2]. 

L2 General Description of the 6-Inch MSBS and ASPS 

Two major magnetic suspension systems are available at Old Dominion University 
(ODU): the 6-inch Magnetic Suspension and Balance System and the Annular Suspension 
and Pointing System. These two systems will now be described briefly. 

In aerodynamic wind tunnel tests, models are normally mounted on mechanical 
struts or stings. Support interference with wind tunnel models continues to be a serious 
problem. For example, it is virtually impossible to make steady-state interference-free 
measurements of wake flow fields without a magnetic support system. MSBS implies 
systems that can both suspend a model (“suspension”) and also measure the forces and 
moments acting on the model (‘balance”). This approach avoids the need for any 
mechanical support of the model. 

The 6-inch MSBS, shown in Figure 1. 1, is a relatively sophisticated and successful 
design. The magnetic field is created by electromagnets located around the test section. 
Measurement of electrical current flowing in each of the electromagnets can be used to 
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determine the forces and moments acting on the suspended model. These forces include 
gravitational, aerodynamic, and inertial loads. The position of the suspended body is 
inherently unstable [3]. Therefore, a closed-loop control system, which includes a position 
sensing system and feedback logic to adjust the magnetic fields, are applied to stabilize the 
position of the body. Further details of the 6-inch MSBS are given in Chapter ID. 

The Annular Suspension and Pointing System, shown in Figure 1.2, is a precision 
payload pointing system designed for the space shuttle. This system consists of several 
sub-systems. The mam sub-system is the Vernier Pointing Assembly (VPA), which utilizes 
magnetic suspension to provide non-contacting isolation and vernier positioning of the 
payload. Five Magnetic Bearing Assemblies (MBA) provide control over three 
translational movements and two transverse rotational mov ement s. Roll control is 
achieved with a segmented A.C. induction motor. Each magnetic bearing assembly is 
equipped with two proximity sensors which measure the position and orientation of the 
rotor. The signals are then fed back to the controller to produce the control signal which 
will be given to the actuators. Further details of the ASPS are given in Chapter V. 

L3 Hardware History and Current Status 

Work on the 6-inch Magnetic Suspension and Balance System was originally 
supported by the Full Scale Research Division of the NASA Langley Research Station 
(now Langley Research Center), Virginia, beginning in 1954, and also by the Aerospace 
Research Laboratories, Wright-Patterson Air Force Base, Ohio, beginning inl957 [4], The 
prototype of the 6-inch MSBS was originally designed by Timothy Stephens [5], Milan 
Vlajinac [6] designed the subsonic wind tunnel which was compatible with the magnetic 
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balance test section in 1966. Construction was completed in 1969 at the Massachusetts 
Institute of Technology (MIT) under contract from NASA Langley Research Center, and 
the unit was used at MIT for research until 1982 [7], The 6-inch MSBS was relocated to 
NASA LaRC in 1984 and resided there for another several years. During these years, the 
magnetic suspension configuration, as originally conceived, did not prove to be completely 
satisfactory in terms of convenience in operation and adaptability to a variety of models 
and test techniques. In September 1989, a program began to restore the system to fully 
operational status with the wind tunnel The basic ideas behind the system remained 
unchanged, but the designs of the various subsystems were refined. The improvements 
included Electromagnetic Position Sensor (EPS) electronics, control system rewiring and 
new water-cooled ballast resistors. Reliability problems still existed, and only limited 
testing was carried out [8]. In 1994, the system was transferred to the Aerospace 
Engineering Department at ODU and is in the first stage of recommissioning and 
refurbishing. 

The history of the ASPS is somewhat similar . The prototype was developed for 
NASA LaRC by Honeywell Satellite Systems in the early 1970’s and delivered to NASA 
LaRC in 1983. However, the system was never assembled at LaRC due to a shift in 
priorities [9], The device was loaned to ODU in late 1992 in a disassembled condition. 
Daniel J. Neff [10, 1 1] carried out the recommissioning. The control electronics rack was 
updated and software was developed to implement A/D and D/A conversions. A 
Proportional-Derivative (PD) controller was written to suspend the Vernier Pointing 
Assembly against gravity in five degrees of freedom 
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L4 Project Goals 

The goal of this project is to continue development of various areas of technology 
in support of the magnetic suspension research program at ODU. The long term focus of 
this research is the recommissioning of the 6-inch MSBS system The first goal involves 
testing the field properties of the MSBS coils and checking against original specifications 
to verify if operation is as designed, and to permit examination of new hardware 
configurations. The next goal is to analyze the design of the EPS and eventually bring it 
back on-line. Finally, a new digital control system must be developed to replace the 
original analog system Since the 6-inch system is not operational, the control system of 
the ASPS has been refined and made somewhat generic through improved software and 
the addition of some hardware. It is expected that this improved controller will later be 
modified and applied to the MSBS. The long term goal is to restore the 6-inch MSBS and 
associated wind tunnel to full operation with a digital controller and a revised magnetic 
configuration. 
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CHAPTER n 
GENERAL SURVEY 


For either a magnetic suspension or a magnetic balance system, an important 
requirement is determination of the magnitude and direction of the unkn own aerodynamic 
forces and torques by solving for them either implicitly or explicitly. In this chapter, some 
basic principles from classical magnetostatics are first presented. The equations for the 
generation of magnetic torques and forces, as well as some special cases, are given. The 
discussion is then extended to the realization of a generalized system. The analysis for a 
specific multi degree of freedom system, the ODU 6-inch MSBS, is given later in 
Chapter IV. 

H.1 Basic Magnetostatics 

A magnetic field H can be generated within a three-dimensional region by some 
external current Un Any ferromagnetic body can be magnetized when immersed in this 
magnetic field. The magnetization M describes the distribution of the magnetic field inside 
the immersed body. The ma gnit ude and direction of M depend upon the shape of the body 
in comparison to the external magnetic field, the history and the strength of the magnetic 
field where the body has been immersed, and the constitutive properties of the material 
from which the body was fabricated. The relationship between the magnetization M and 
the total magnetic flux density B (B = pH) of the field around the body may be described 
as [12]: 

(D-l) 


B = poM-Kp-po )[D]M = (po[I]+(p-Po )[D])M 
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n-wd+x) (n-2) 

where [D] is the demagnetization factor determined by the body geometry, Mo is the 
permeability of free space, p is the permeability of the actual medium (assuming 
M«r = Mo), and x is the magnetic susceptibility. 

Different materials have different abilities to generate internal magnetic fields. 
Diamagnetic materials have negative susceptibility and appear to exclude the external 
magnetic field, while paramagnetic materials have positive susceptibility and appear to 
“attract” the field. For ferromagnetics, there is a highly non-linear relationship between the 
induced magnetic field and the applied magnetic field, and the induced field has hysteresis 
in the sense that the field is non-zero after the applied field is reduced to zero. 

For a typical ferromagnetic body m»Mo and |xD;p» 1, therefore Equation (II- 1) 
can be written as: 

[D]M » H (H-3) 

This approximation implies that the magnetization level in the body is essentially 
independent of p in the linear range of magnetization. Therefore, M is only determined by 
the material and shape of the body placed inside the field H where 

M = [D] ‘H 



"m x " 


^XX ^XY ^XZ 

-I 

H x " 

or 

m y 

= 

DyX Dyy DyZ 


h y 


M z 


_Dzx Dzy Dzz 


1 

X 

N 


In the application to wind tunnel systems, two sets of axes are in common use. 

One is the balance axes X, Y, and Z, the other is the principal magnetic axes a, b, and c of 
a non-axisymmetric body. These two sets of axes are presented in Figure II. 1. Notice 
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these two sets of axes are related by Euler angles of \j / (yaw angle), 0 (pitch angle) and cp 
(roll angle). A transformation matrix [R] ([a b c] T = [R][X Y Z] T ) can be defined by these 
three angles. The development of [R] is given in Appendix A. 

In the body coordinates a, b, and c, Equation (0-4) becomes: 


X 


D. 

o -L 

D b 


X' 

M b 

= 

0 

H b 

. M c. 



1 

H 0 . 


0 




(n-5) 


The demagnetizing fiict ors are related to each other as follows [12]: 

D.+Db+Dc = 1 

By using wind tunnel coordinates. Equation (II-4) becomes: 


1 

5 

X 

1 


1 

X 

X 

1 

m y 


h y 

M z 


1 

X 

N 


(H-6) 


(n-7) 


H, 


[[»]■'[!>] m] |Hv 

H, 


-nr 


D. 

0 
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D fc 


H, 


[nrnw'n]] h y 

H, 


D. 



(n-8) 
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IL2 Force and Torque 

When immersed in a uniform magnetic field, an unconstrained magnetized body 
will experience a torque which turns the body to align its magnetization axis with the 
direction of the applied field. This behavior is known as the compass needle phenomenon. 
In reverse, a torque can be generated by rotating the body away from its equilibrium 
position at a certain angle. The torque can be represented as [12]: 

T = VMxB (0-9) 

for a uniform and steady flux density B, or 

T = J M x B dV ( n - 10 ) 

for a non-uniform steady flux density B. Here V is the volume of the levitated body. 

Some special cases are discussed when only uniform B is involved [5], 
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M x B = 


i 

A 

♦ 

J 

k 

H. 

H b 

H„ 

D. 

D, 


mh. 


mh. 


H,B. H.B, 


D, 

H.B b 


D. 

H b B. 


D. 


D. 


(n-i4) 


Case c. The body has equal demagnetizing factors (a spherical shape) and hysteresis 
effects are negligible. 

D, = Db= D c = 1/3 = D 
H. 


M = [D]"‘H = 


D 

ILl 

D 

JL 

D 


= — H 

D 


(D-15) 


MxB = — HxH = 0 
D 


(n-i6) 


Case d. Applied field is parallel to any of the principal magnetic axes. H = H, i and 

M = [D]‘H= HJD.i 


M x B = 


1 

H. 


j k 

0 0 


D. 

pH, 0 0 


(II- 17) 


= 0 


In cases a and b only pitching and yawing moment can be generated (T, = 0). In cases c 
and d, torque can not be generated along any axis (T — 0). 

The generation of rolling moment has always been difficult in wind tunnel magnetic 
suspension systems. Stephens [12, 13] suggests three concepts of roll moment generation: 
a) Interaction between a transverse or non- symmetrical model magnetization and 
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transverse D.C. field components, b) Interaction between a controlled single phase 
transverse A.C. field with a model supporting a conduction plane or loop, which perhaps 
can be located inside the wing area, and c) Interaction between a rotation 2-phase 
transverse A.C. field and a copper plated or other conducting model. 

The only way to generate a net force on a magnetized body in steady fields is by 
applying a non-uniform B. In a Cartesian coordinate system [12], 

F = J (M • V)B dV ( f o r magnetized body) (II- 1 8) 

F = VM • VB (for permanent magnet) (11-19) 


Notice that because VxB = 0 (Maxwell’s Equations), 


F = 


5B V 
M x — - 
x dx 

3B V 


dB v 
+ M x 
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0X 
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dx 
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d z 


dB Y 
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dx 

dBy 
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+ M, 
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dy 
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d z 


3B y 
Y dx 

+My ^ 

ay 

v dz 


+ M Z ^ 

z dz 

d b, 

+M ‘ar. 

3B Z ' 

3x 
dB z 

dy 
dB z 


+ M ; 


+ M. 


+ M. 


(n-20) 


dz J 


If magnetization only exists along the X direction, which is also the body’s axis of 
symmetry, the force equation can be simplified and forces can be controlled independently 
by adjustment of the gradient of B: 
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F = 


3B v 
M x — 2L 

x dx 

5B V 

Mx ~dy" 

M X ^- 

dz 


(n-2i) 


H3 System Analysis 

In general, forces and moments will be generated on the model Therefore, in a real 
system, application of the above concepts relies upon the ability to independently control 
and compensate the individual force and torque components. A successful system not only 
includes features such as a source of magnetization and a source of field gradients, but 
also needs several subsystems to complete the suspension task. These subsystems include 
power supplies, a cooling system, and a sensing and control system to maintain the 
stability of the suspended object and overcome couplings among multiple degrees of 
freedom. The more degrees of freedom the system has, the more complicated the problem 
becomes. In this section, force and torque generation, and the corresponding control law 
for a one degree of freedom system analysis, is given first as a microcosm of a complete 
system. General descriptions of sensing units, power supplies and cooling systems for 
multi degree of freedom systems are also given. 

EL3.1 Force and Torque 

The one degree of freedom system to be discussed consists of a single N-tum 
si mila r coil of radius a which is shown in Figure II. 2. The magnetic field on axis is given 
by 


H z = 


Nla 


2 * 3/2 


4rt(z + a ) 




Nla 2 


H, 


2(z 2 + a 2 ) 3/I (z 2 / a 2 + 1) 


3/2 


(0-22) 
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where Ho=NI/2a 


dH z _ 3z H 0 (n-23) 

dz a 2 (1 + z 2 /a 2 ) 5/2 


F z=K^ dV = /MM 2 ^dV 


dz 

r 3pz M Z H 0 
' a 2 (1 + z 2 / a 2 ) 


dV 


(n-24) 


If Mz is generated by a permanent magnet, it can be assumed to be a constant. Therefore 

3[iM z V zH 0 


F z =- 


a 2 (1 + z 2 /a 2 )” 2 


(0-25) 


is zero at z = 0 and is ma ximum at z = ^ a . If Mz is generated by magnetizing a 


ferromagnetic material in the linear range, then Mz = Hz/D, where D = D, = Db = D c = 
1/3. Hence, the force component is given by: 


F z = JM* VB z dV = J|i — ^-dV 


2/3p,z 


= -(H 0 ) 2 (^) 


D dz 
V 


aD' a (1 + z 2 /a 2 ) 4 


is zero at z = 0 and is maximum at z = — a 

7 


(0-26) 


0.3.2 Control Law 

To study the possibility of achieving a stable system, it is necessary to study the 
motion of the sphere from Newton’s second law: 


z = g + F M /ni 

Zo is defined as the equilibrium position when v = z = 0, z = 0 and mg = -F 0 


(H-27) 
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The control law requires that F increase when the sphere falls and vice versa. A 
Proportional-Derivative controller is defined by: 

H = H(z,v) = H 0 + K,Z + K 2 Z (H-28) 

Notice that in Equation (11-28), Z (= z-zb) is the displacement from the equilibrium point 
and Z is the time derivative of Z which is equal to the velocity. 

mZ = -F z + F m = AF (n-29) 


Substituting (11-25) into the Equation (11-29), 


dF = (- 
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and hence 
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Therefore, the equation of motion for the system can be written as 
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The system dynamics are uniquely determined by the eigenvalues of the matrix on 
the right-hand side of Equation (11-30). 

7 . L,K 2 ^ K,L, — L 2 


Z + -^-Z+- 

m m 


-Z=0 


(n-31) 


The natural frequency is then 


co. = 


K,Lj - L 2 


m 


(n-32) 


and the damping factor is 


c = 


K 2 L t 


2 V m(K,L, - Lj) 


(H-33) 


The control law given by Equation (11-28) is adequate for a pure suspension 
system, but not for the case where a drag load is applied. With a drag load. Equation (II- 
31) becomes: 

a.-IAt+ b .r -Krfx z+P < n - 34 > 


m 


m m 


where Z '= at Z = Z = 0. 

0 L 2 - K,L, 

This new equilibrium position is determined by the magnitude of drag, and the 
steady state error caused by this additional force cannot be compensated for by using a PD 
controller. To eliminate this unfavorable displacement and hold the model at a fixed 
position, the control law is modified by adding a new state variable u which is defined as 
u = j Zdt to get a Proportional- Integral-Derivative controller [12] such that 
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EL3.3 Position Sensor 

To use an automatic control system, position sensing equipment is necessary. The 

desired characteristics for a position sensing system are [12]: 

(a) No mechanical contact with the body exists. 

(b) Translational displacements of the model parallel to the orthogonal wind tunnel axes 
will be transformed into three independent high-level electrical signals proportional to 
the translational displacements over a suitable range. 

(c) Angular displacements of the model about the wind tunnel axes will be transformed 
into three or more independent, high-level electrical signals proportional to or strong 
functions of the angular displacements, over a suitable range. 

(d) Resolution, repeatability, and frequency response must be compatible with desired 
accuracy of static and dynamic wind tunnel results. 

(e) The system must be easily adaptable to a broad range of model geometries. 

(f) The system must be compact, and must not interfere with access to the test section or 
prevent visual contact with the model 

(g) The system must operate without interference from the magnetic fields used to support 
the model. 

Two methods have been applied to detect the model positions. One is optical 

sensing, another is electromagnetic position sensing. Optical systems meet most of the 
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above requirements and are dominant at the present time. A basic system, shown in Figure 
0.3, is needed for each degree of freedom being controlled. The large beam system [12] 
permits the model to be controlled over a large distance. This system is shown in Figure 
0.4 The use of triangular mirrors can double the sensitivity of the system 

The Electromagnetic Position Sensor (EPS) is believed to have low sensitivity due 
to the large air gap and relatively small model size. However, a well-designed, compact 
EPS system, of which only one exists so far, has been demonstrated. The advantage of the 
EPS is its capability to accommodate a wide range of geometry with minor adjustment. 

The system is based upon differential transformer action. The idea is that in an AC. 
magnetic field, the model’s displacement changes the coupling between excitation 
windings and secondary coils. This signal can be picked up, demodulated and decoupled 
to resolve the model’s movement. The details of the hardware arrangement and analysis of 
this system are discussed in Chapter HI. 
n.2.4 Power and Cooling 

The power supplies are usually the most expensive component of the entire 
magnetic suspension and b alan ce system Also, many of the reliability problems and 
performance limita tions are directly related to the power supply system A fast response is 
required which means the power unit must have the ability to absorb the regenerated 
energy when inverting the output. Safety and reliability are also considered to be of prime 
importance. If using an EPS, there is an additional task of minimizing high-frequency noise 
to prevent interferences between two magnetic fields. Two passive filter configurations are 
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shown in Figure 0.5. These low-pass or band-pass filters can be connected after the power 
supplies to limit the output frequency. 

The suspension coils must be designed to handle large ohmic power dissipation. 
Three schemes have been used for heat removal purposes. They are an air-cooled system 
with fans added, a conduction cooling system (shown in Figure n.6), and a liquid cooling 
system. The liquid used in a liquid cooling system could be either water or oil. The liquid 
usually runs in a closed loop to minimiz e the requirement for filtering, dechlorination and 
other liquid treatment. If water is used, it must be demineralized and deaerated. Some 
attached subsystems such as pressure, temperature and flow indicators are needed to 
protect the system from overheating. A complete description of a water cooling system 
used in the 6-inch MSBS is introduced in Chapter HI. 



Rotation Sequence: X 

vy - about Z 
9 - about new Y 
<t> - about new X 
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Figure IL1 Model and Wind Tunnel Axis System 



N turns 


Figure 11.2 One Degree of Freedom System 
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Figure 0.4 Large Beam Position Sensor 
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CHAPTER ID 

HARDWARE DESCRIPTION OF THE 6-INCH MSBS 


A detailed description of the 6-inch Magnetic Suspension and Balance System is 
given in this chapter. The configuration of the electromagnetic coils, which generate 
controlled magnetic fields in the test section, is introduced first. The general features of 
the model position sensing system, compensation system, and power and cooling systems 
are then outlined. Current status of the 6-inch MSBS system is presented at the end of this 
chapter. 

HL1 Magnet System 

Although a variety of magnet configurations can be conceived, the design should 
be as compact as possible for reasons of power economy. The geometric constraints 
imposed by the wind tunnel, the position of the EPS, the accessibility to the test section, 
and the view path for a Schheren system must also be considered. Further requirements 
include providing access for installation of in strumentation probes and injection of the 
model. 

The prototype magnet system hardware is comprised of several elements and 
subelements. The magnetic fields are generated by the Helmholtz coils, saddle coils and 
iron core magnet assemblies. Other structures such as end plates, tie rods, and safety 
shrouds are alcn necessary to support the coils, to provide the structure, and to protect 
other devices. These subsystems are now introduced in turn. 
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m.1.1 Helmholtz Coils 

Helmholtz coils provide an axial magnetization field and drag force control. A 
“Helmholtz configuration” is defined as a pair of identical circular solenoid coils between 
which the axial spacing is approximately one mean coil radius. The Helmholtz coil system 
used in the 6-inch MSBS has two subelements: magnetization windings to provide a 
steady state magnetic field, B x , and drag windings to provide a gradient of the axial field, 
Bxx- 

Figure III. la shows the arrangement of the Helmholtz coil system. Each coil is 
wound in three parts which are separated by insulating layers. The middle sections of the 
two coils are arranged in series to generate uniform B x over an appreciable volume of the 
central region of the test section. The inner and outer sections are connected in series, 
with the current in the two coils circulating in the opposite sense. Therefore, at the 
centroid of the test section, the axial field strength is zero, and a uniform gradient Bxx is 
produced along the axis. 

The coils are wound of hollow copper conductors in layers of constant radius. 
Joints between layers occur at the be ginning and end of each layer, and they incorporate 
hydraulic fittings which allow cooling water to be circulated. Figure HI. lb shows the coils 
before installation [5]. 

HL1.2 Saddle Coils 

The transverse magn etic field components By and Bz are produced by two pairs of 
saddle shaped coils. The arrangement and a photo are shown in Figure HI. 2a and HI. 2b. 
These coils consist of an inner pair and an outer pair which are orthogonal to each other. 
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By and Bz are then composed of components from both pairs of coils to provide pitch and 
yaw moments. The apertures in the inner coils overlap the apertures in the outer coils in 
order to provide a horizontal viewing and Schlieren path, and a vertical path for additional 
access to the test section is available. 

The coils are wound of hollow copper conductor in double layer pancakes, bent to 
shape and encapsulated. Cooling manifolds are attached to the ends of coils [5]. 

HL1.3 Lift and Side Force Coil Assemblies 

Control of lift and side force is affected by a pair of iron core magnet assemblies, 
shown in Figure m.3a. Each assembly consists of octagonal laminated silicon steel core 
structures having four symmetric tapered poles which point toward each other and are 
canted at 45° to the plane of the octagon. The iron core is used to enhance the magnetic 
field strength by an appreciable factor thereby reducing the power consumption, which is a 
highly desirable characteristic. On the other hand, the iron places a limit on the range of 
linear operation of the magnet system, due to saturation effects. Furthermore, the 
inductance of the coil system wifi increase due to the presence of the iron cores. This 
increment will generate more reactive power in the circuit and nearly offset the power 
saved by the merit of strength enhancement [4]. 

Four identical magnet coils are mounted on the joints of the poles vertically and 
horizontally. The coils are constructed of hollow copper conductor, which is water- 
cooled. The four coils with vertical axes (two coils per assembly) are connected in series 
and control the lift force, shown in Figure ID. 3b. The four coils with horizontal axes (two 
coils per assembly) are connected in series to control side force. A picture of these 
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assemblies before installation is shown in Figure m.3c. 
in. 1.4 Roll Control System 

A complete roll control system compatible with the 6-inch MSBS system was 
originally suggested by Stephens [5] to provide control of the sixth degree of freedom, 
roll. This system will not be introduced here since it relies on A.C. fields to provide 
separation from the other five degrees of freedom The development of this system, which 
includes the drive circuit, the Roll Position Sensor (RPS), and the compensation circuit, 
was completed in 1982 [14, 15]. 

EDL1.5 Mechanical Considerations 

All magnet coils are provided with thermal sensor switches and coolant flow 
switches to prevent thermal overload. The various coils in the magnet system are 
interconnected and connect to terminal boxes on top of the assembly. The magnet 
assembly is sandwiched between two fabricated cruciform aluminum end plates which are 
held together by longitudinal tie rods. Lateral alignment of the subsystem is provided by 
inserting dowels into structural elements. Figure m.4a shows the structure of the magnet 
system, and Figure m.4b shows a photo of the structure before assembly. The complete 
system is then covered by safety shrouds which prevent personnel from contacting high 
voltages in the magnet coils and to prevent magnetic objects from being attracted into the 
magnet system which might cause short-circuits. A support structure is also designed to 
distribute the weight of the system uniformly to the floor. Finally, the entire system is 
connected to the 6-inch wind tunnel with an octagonal test section, shown in Figure m.5. 
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DL1.6 Wind Tunnel 

The wind tunnel compatible with the magnetic balance system is an open circuit, 
closed-jet tunnel The contraction is constructed of molded plywood. Two settling 
chambers and a honeycomb are located upstream of the contraction to improve flow 
quality. The octagonal test section is located within the magnetic assembly. It is made of 
transparent acrylic, measures 6.25 inches across flats and is 36 inches long. A two-piece 
diffuser with access doors is located between the magnetic assembly and the squirrel cage 
fan. A ma ximum dynamic pressure of 2.5 psi, a Mach number of 0.6, and a freestream unit 
Reynolds number of 3.5xl0 6 /ft are attainable [6]. 

HL2 Electromagnetic Position Sensor 

Although optical model sensing systems have been used extensively because of 
their linearity and low noise, the most obvious drawback of this approach is that the 
design of an optical system must be tailored to the size and configuration of the wind 
tunnel modeL The Electromagnetic Position Sensor is an alternative, which employs the 
principle of the differential transformer with the model forming the core to resolve five 
components of the model position and orientation. 

An EPS compatible with the 6-inch magnet system was originally designed by 
Stephens [16]. The simplified schematic of the EPS system is shown in Figure HI. 6. The 
sensor subsystem, which is a transducer coil assembly, and the electronics subsystem, 
which is used to actuate the coil system and resolve the output signal, are introduced 
respectively. A qualitative description of the system's operation is also given. 
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m.2.1 EPS Coil System 

The key part of the overall system, (e.g. a coil assembly with excitation and pickup 
coils) is shown in more detail in Figures m.7a, HI. 7b and m.7c. The excitation and pickup 
coils are arranged to provide a viewing window into the test section. The windings are 
potted with epoxy for mechanical strength and painted black outside and white inside to 
enhance the model visibility during operation. 

The excitation windings consist of a pair of Helmholtz coils, designated as EXC A 
and EXC B, which terminate at the end of the assembly at a 2-pin socket. These coils 
generate a uniform axial magnetic field which oscillates at f = 20 kHz when driven by a 
suitable power amplifier. 

The pickup coils consist of seven pairs of No. 30 Teflon-covered stranded small 
gage wire, all having an equal number of turns (thirty turns on each coil) and terminating 
at a multiple-pin socket at the end of the assembly. Axial pickup coils, designated as AX A 
and AX B, are also in a Helmholtz configuration and are coaxial with the excitation coils 
but are connected in opposition. These coils are used to pick up the model’s axial signal. 
The remaining six pairs of pickup coils are in four rows of saddle-shaped coils, from which 
signals are resolved to derive the model’s horizontal and vertical translation, as well as 
pitch and yaw orientation. The parameters of the pickup coils are given in Table HI. I [15] 
and compared with the data obtained from the pickup socket. 

A small electrical box, connected directly behind the excitation signal cable, is used 
to minimize the reactive power. The circuit is shown in Figure EH 8. The required 
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Table HL1 ODU 6-Inch MSBS, EPS Coil Parameters 



Rrf= 6.2 n 
R*= 11.7 £2 
Rd = 14.0 £2 
R*. =9.5 £2 
R«= 8.5 £2 

Rji = 6.6 £2 

R ^ = 6.2 £2 


Pickup Socket 


EPS Coil Circuit 

DC Resistance 1 

m 

Inductance 2 

(mH) 

Self-capacitances 3 

(mF) 

Axial 

14.76 

0.919 

0.00141 

Lateral I 

7.57 

0.468 

0.00156 

Lateral II 

10.41 

0.517 

0.00143 

Lateral HI 

7.18 

0.416 

0.00126 

Vertical I 

12.42 

0.509 

0.00144 

Vertical II 

9.49 

0.502 

0.00146 

Vertical in 

12.47 

0.488 

0.00152 


1. Measured with Fluke digital VOM. 

2. C alculat ed from resonance frequency with 0.25 pf mica transmitting capacitor. 

3. Calculated from self-resonant frequency. 





30 


capacitance to offset the inductance of excitation coils can be estimated as follows: 


— (sL + R) 




sL + R 


(sL+R) + — s 2 LC + RCs+ 1 
Cs 


(m-i) 


S —►JO 


(R + jcoLXl - LC© 2 - RC©j) 

''Opmioaimptdnoa " (J _ LC( o J ) Z + (RC©) 2 

_ R + j(©L-L 2 C© 3 -R 2 C©) 
(1- LC© 2 ) 2 + (RC©) 2 


(ffl-2) 


To minimiz e reactive power. 


©L-L 2 C© J -R 2 C© ->0 


Therefore, 

c k 

o 2 L 2 +R 2 (DI-3) 

For L = 404 (iH, R = 0.7 £2, and o = 2nf= 4rcx 10 4 , Equation (HI- 3) gives 
C = 0. 1567x 10 -6 F. The capacitance actually used is 0.2 xlO -6 F. 

The working principle of the coil assembly is briefly introduced here. Once the 
excitation coils are actuated, they produce a uniform magnetic field over a region 
surrounding the centroid of the test section, oscillating axially at 20 kHz. When a slender 
ferromagnetic model is aligned with the axis of the coil assembly, it will be magnetized at 
the same frequency with the excitation field. Figure IH.9 shows the induced fields when 
the model is at different locations in the frame. From the diagrams, it follows that the drag 
signals are measured by AX without first-order coupling with the other four motions. A 
general statement can be made that linear motions are detected by subtracting coil outputs 
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and rotations are detected by summing coil outputs. For example, vertical motion can be 
resolved by subtracting the output of pick-up coils V3 from the output of VI, and pitch 
motion can be resolved by s ummin g VI, V2 and V3 together. The results apply equally 
well to lateral translation and yaw rotation by symmetry. 

IIL2.2 Electronics System 

Figure m. 10 shows the block diagram of the electronics associated with the EPS 
coil system. However, this original system provided poor sensitivity, high noise and 
correlated interference, some of which is due to the weak points in the signal processing 
circuit design, and some due to the unreliable analog elements. 

A new design was developed in the late 1980s and the old data reduction circuit 
was replaced by a new compact electronic system providing output with less noise, drift, 
and higher sensitivity and repeatability [17], EPS circuit diagrams are given in Appendix 
B. Unfortunately, system documentation was never completed. Therefore, some efforts to 
analyze the circuits were made and the comments are given in Appendix B. The 
compensation system consists of several stages in cascade, and each of them has several 
elements in co mm on Figure ED. 1 1 is a block diagram of one channel, and it will be used 
to illustrate the operation of the whole system. Other channels are similar. 

EDL2.2.1 EPS Oscillator and Power Amplifier 

The 20 kHz, 2.5 V peak-to-peak carrier sinewave for primary coil excitation, input 
stage n ullin g and demodulation is supplied by a crystal-controlled IC oscillator with 
±0.05% stability and < 5% total harmonic distortion. The original vacuum tube amplifier 
was replaced by a new solid state power amplifier. This item did not arrive at ODU with 



other EPS hardware. According to the documentation [17], this amplifier should be an 
audio frequency amplifier with output amplitude about 210 vohs peak-to-peak. 
m.2.2.2 Reference Phase Adjustment 

The output of the amplifier is sent to both the position sensor coil assembly and the 
nulling system as the null-reference used in the demodulation stage. Signals pass through 
amplifiers driving bipolar balanced output transformers to provide seven pairs of active 
summing networks and five demodulator references. 

HL2.2.3 Summation and Demodulation 

The seven differential coil voltages are sent to band-pass filters and amplified 
before arriving at s ummin g circuits, in order to attenuate higher order harmonics. Five 
AD630 demodulator integrated circuits combined with low-pass filters are applied to 
obtain continuous DC outputs which should be approximately proportional to the 
movements of the model relative to the axes of the transducer assembly. 

DL3 Compensation System 

Since the suspension system is a position control servo, designed to maintain the 
suspended model at some co mman ded position, a superior compensator should 
accomplish at least three things: 1) Stabilization of the loop, 2) Minimization of static 
position errors for given disturbance inputs, and 3) Minimization of static position errors 
for given position co mman d inputs. 

A functio nal block diagram of the suspension system is shown in Figure III. 12. The 
position of the model is translated into electrical signals by means of the position sensing 
device, and these signals are compared with the position command signals. The resulting 
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position error signals are modified by passing through a set of feedback compensation 
networks, and then amplified and applied to the actuators. The magnetic forces and 
moments, which are functions of the coil currents, act on the model to cause changes to its 
position and orientation. These operations adjust continuously to counteract the 
aerodynamic and gravity loads. 

The feasibility of applying digital information processing to the magnetic 
suspension system was proposed in the mid-1970s [18, 19]. Establishment of the overall 
system model and calibration procedures for discrete subsystems were accomplished 
during this period. Although a more advanced controller might increase operating 
performance levels for this non-linear control problem, the control method outlined in 
References [18, 19] provides a framework for further recommissioning work. 

A simple and effective form of the controller consists of three paths: proportional, 
integral and third-order lead-lag network, shown in Figure m. 13. The main purpose of 
using the integrator is to improve the system’s low frequency response, in other words, 
maintain zero position error with a steady disturbance input. The application of a 
proportional controller helps reduce rise time, and the lead-lag path helps reduce 
overshoot. From a practical point of view, experiments are necessary to determine 
acceptable combinations of settings for the gains and poles and zeros of the compensator. 
These settings are dependent on the model used. 
m.4 Power and Cooling System 

The power system is the most expensive part of the 6-inch MSBS. This subsystem 
was modified several times to increase drive current capability [4, 5, and 20]. A schematic 
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of this old power supply system is given in Figure III. 14. Seven power amplifiers are 
needed for this six degree of freedom system, Le., five D.C. amplifiers for drag, lift, side 
force, pitch and yaw respectively, and two additional A.C. amplifiers for rolling moment 
control. The D.C. amplifiers are required to operate in all four quadrants. Therefore, they 
could provide a continuously variable current output which is linearly proportional to the 
input voltage control signal. The repeatability and stability of the output current is of the 
order of 0. 1% of the operating output current for day-to-day operation. The maximum 
output power level is approximately 1,000 amps into a 0.34 (±20%) ohm load. The L/R 
time-constant of this load was about one second. The power at maximum resistance is 
approximately 400 kW. 

The magnet coils, and the conductors which form the interconnections between 
coils, must be provided with a means of removing the heat which is generated due to 
ohmic losses. The cooling system must be reliable and of moderately high power-density 
capability. The most suitable choice of cooling media appears to be water, and in order to 
avoid electrolysis and silting of the conductors, the cooling water needs to be 
demineralized and deaerated. 

Figure HL 15 gives the cooling as well as electrical paths in the conductor coils. 
The brass adapters fit nylon tubing (d = 0.25 inches), which is used to connect the coils to 
the inlet and outlet coolant manifolds. The cooling water actually circulates through the 
magnet coils and is contained in a closed loop, as shown in Figure HI. 16. Heat is removed 
from this water with a heat exchanger, which transfers the heat to unpurified water. This 
unpurified water is then circulated through a cooling tower. The “inner” cooling loop 
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(purified water loop) is provided with a high-pressure circulating pump to force the 
purified water through the magnet coils. A reservoir is used to accommodate changes in 
volume of the coolant due to expansion, and a portion of the cooling fluid flow is 
circulated through a small demineralizer and deaerator to main tain the required purity. 
Notice that pressure, temperature and flow indicators are also necessary to prevent 
coolant failure which will result in magnet coil burnout. 
m.5 Current Status of the 6-Inch MSBS 

The 6-inch Magnetic Suspension and Balance System was transferred to the 
Aerospace Engineering Department of Old Dominion University in 1994 and is in the first 
stage of recommissioning work. 

The magnet system was checked and tested to verify the magnetic field 
distribution. Results are given and discussed in the next chapter. A new amplifier is needed 
to complete the EPS hardware. Further, schematics of the printed circuit boards, which 
are very helpful to recover the EPS electronics, are sought, but have not been found at this 
time The original power supply and compensation systems have been totally abandoned 
due to their out-of-date condition. These two systems need to be redesigned or replaced 
according to new performance requirements. New technology in digital electronics should 
bring a more compact and advanced control system. The aged cooling system needs to be 
either repaired or replaced. Restoring the whole wind tunnel is the last task of the 
reassembly work. 
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Figure HL2a Saddle Coils: System Arrangement 



Figure IIL2b Saddle Coils: Prior to Installation 
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Figure IIL3a Lift and Side Force Coils: System Arrangement 
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Figure IIL3b Lift and Side Force Coils: Flux Paths in Lift Force System 
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Figure DL3c Lift and Side Force Coils: Prior to Installation 
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Figure EDL5 A Schematic of the 6-Inch MSBS Wind Tunnel 
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Figure ED.6 A Simplified Schematic of the EPS Sensor 
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Figure HL7c EPS Coil Assembly; Prior to Inst3.lia.tion 
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Figure IH.8 Circuit Used to Minimize the Reactive Power in EPS Coil Assembly 
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Figure III.9 EPS Signal Contributions from Single Flux Lines 
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Figure m.10 EPS Electronics System Block Diagram 
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Figure III. 11 Block Diagram Comparison of Pitch Channel 
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Figure in. 12 Functional Block Diagram of System 



Figure HI. 13 Schematic of Compensator 
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Figure HI. 14 Power Supply System Schematic 



Figure 1H.15 Layer-wound Hollow Conductor Coil Cooling and Electrical Paths 
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CHAPTER IV 

THE MAGNETIC FIELD OF THE 6-INCH MSBS 


In this chapter, the magnetic field generated by the magnetic assembly of the 
6-inch MSBS is estimated by analyzing the individual subsystems. Measurements and 
computational results are also given and compared with each other. 

IV. 1 Performance Estimation 

In the design of a large-scale magnet system, it is desirable that the magnetic 
performance be estimated accurately. Two classical approaches for estimation are 
commonly used: 1) the use of a small scale model and 2) the analytical calculation of the 
magnetic fields. The former is more accurate and more reliable for a complicated system, 
but it also more expensive. The latter is much cheaper, although the simplification of the 
systems’ geometry, saturation effects and the absence of the iron cores may cause large 
errors. Modem computational techniques are available and are superior to both of these 
approaches, but require the use of sophisticated computer codes. 

The basic formula used to calculate the magnetic field generated by coils with 
different geometries is the Biot-Savart Law, which states that the magnetic field B at a 
distance r from a line current element Idl is: 


B 


= J >-1 


pldl r 

r * — 

27tr r 


(IV-1) 


IV. 1.1 Helmholtz Coils 

The calculation begins with a single loop shown in Figure IV. 1. The magnetic field 


generated along the Z axis is Bz. 
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= Ugla 2 ft. = Mpla 2 (IV-2) 

z 4 tc(z 2 + a 2 ) 3/2 J o 9 2(z 2 + a 2 ) J/2 

For the Helmholtz coils in Figure IV.2, the magnetic field at the center of the axis 
of symmetry is the superposition of the fields generated by the drag windings and 
magnetizing windings. The fore-and-aft pair is used to improve field uniformity. 

The symmetrical part Bx , generated by the middle sections of the Helmholtz coils, 
is the main axial magnetization, called the “bias ” field. The current per unit area in the 
windings is 

K„ = — 

(x, -x t Xr, -r 2 ) 


B x (0,0,0) = (0,0,0) + 8^ (0,0,0) = 2B 

Xkft (0,0,0) 

Ho r ’K 0 ^«r 2 K 


= ff M N dS= pp .y ^ dxdr 

" (x +r ) V«i (x +r ) 

r , +V r 3 2 +x i 2 


= HoK, 


r 


r 2 + x 2 2 


x 2 In 


r, +^/fj 2 + x 2 


(x 2 - XjXtj - r 2 ) v. 2 r 2 +^/r 2 2 +x 2 


= - x. In 
1 * 




+ X 


(IV-3) 


i J 


where po = 4 tix 10' 7 henrys/meter, r 3 = 14 inches, x 2 = 8 inches 

N = 400, r 2 =10 inches, xi = 4 inches 

Substituting these numerical values into Equation (IV-3) and converting to meters 
(multiply by k — 39.37 inches/meter) to get the correlation between current and magnetic 


field yields 
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Bx/I — 1 1.78 gausses/amp = 1. 178x 10’ 3 tesla s/amp 


(IV-4) 


5B 

The drag force is produced by the axial gradient field Bxx ( = — — ) which is the 

dx 

anti- symmetric contribution. Defining 

NI = Ko(r 4 -r 3 ) = K 0 ' (r^rO 
Xo=(Xi+X2)/2 


n _ r f 4 j , ff 2 r j 

B ** ~ K 7Z — 72 , 3/2 ' + K 2 ~ ~2 ) 3/2 * 


( x o +r ) 


(IV-5) 


B _ p_3M,K 0 r‘x ldf+ p_ M ,K/r i x, di . 
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= -3n.x.[(-^-X- _ ttA +< N,t 


r 4 -r 5 3x/(x 0 '+r 2 r* ir, 

3 


•X — 

2 /„ 2 


o v^o ' * / i*4 r 2 r i 3x 0 (x 0 + r ) 


2 \3/2 


r / Njl y r* 13 x 

- -^x„[(— X x p (x<1 + r p )M - x p (x p +rjl) » ) 

+(— 2 — -X 


r 2 - r i x 0 2 (x 0 2 + r 2 2 ) 3/2 x 0 2 (x/ + r," ) 3 


)l 

2,2 v3/2 ^ 


(IV-6) 


where Ni = 200, r L = 8 inches, 

N 2 = 200, r 2 = 10 inches. 


r 3 = 14 inches, 
r 4 = 16 inches. 


xi = 4 inches 
X 2 = 8 inches 


Bxx /I = 1.34 gausses/inch/amp = 5.27x10 3 tesla s/m/amp 


(IV-7) 


IV. 1.2 Saddle Coils 


Using the Biot-Savart Law for the geometry shown in Figure IV.3, the magnetic 
field at point P is 

B = 7 — T C 0 S P d P e e = 7 ^(sinp 2 -sm3,)e e . 

4xa J Pi 4xa 


(IV-8) 
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For the geometry shown in Figure IV. 4, the magnetic field generated by a current element 
at point O is 

B z = -J "’* 2 — — sin <|> sin pd<J> 

z 4u(a 2 + L 2 ) Y 

= -^-cos 2 Psin^cos^j +cos<|> 1 ) (IV-9) 

47ta 

The field strength of the saddle coil is complicated to calculate due to its geometry. 
However, an approximation can be given by using a “mean turn”, corresponding to the 
centroid of the winding, as shown in Figure IV. 5. The transverse field B t , at the origin, can 
been seen as the contribution from the straight line currents B t i and the contribution from 
the current in the segment of the loops B a . 

Applying Equations (IV-8) and (IV-9) to calculate B t i and Ba, and noting that 

L 

<|>i = <Ji 2 , sin 0j = - sin P j = — - — , and that each coil consists of two saddles, 

K~) J +r 3 ] 1/2 


yields 


B. = 


pNI 

4kt 


(1 + 


[1 + 


(L/2) 


j] 


VI 


l+(— ) 2 


) x 4 x 2 cos<j> 


For the inner saddle coil, 

B^in)/I = 5.5 gausses/ amp = 0.00055 teslas/amp 
where N = 88, r = 5.2 inches, L/2 = 5.5 inches, (J> = 15° 

For the outer saddle coil, 

Bt(o U t)/I = 5.58 gausses/ amp = 0.000558 teslas/amp 


(IV- 10) 


(IV- 11) 


(IV- 12) 


where N = 133, 


r = 7 inches. 


L/2 = 5 inches. 


4> = 15° 
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IV.1.3 Lift and Side Force Coil Assemblies 

A lift coil model is shown in Figure IV. 6 to calculate the magnetic field strength in 

the Z direction along the wind tunnel axis. Using the Biot-Savart Law, the field is 

B _ (iNI |*J" dl x s 
271 J 0 s 3 

Noting that 

s = r-a 

= ri-acos0i-asin0j 


s 3 = [(r - a cos0) 2 + (a sin0) 2 ] 3/2 
dl = ad0(- an0i + cos0j) 


dl x s = a(a - r cos0)k 


the Z component of the field is 


_ pNI r 2 » a(a - r cos0)d0 

2 2i t (r 2 +a 2 -2arcos0) 3/2 (IV- 13) 

where N = 290, a = 3 inches, r = 11.5 inches 

For a pair of coils, arranged as shown in Figure IV. 7a, Bz along the X axis of the 
tunnel is calculated and shown in Figure IV. 7b. The slope of the magnetic field per amp is: 


Bxz /I = 0. 13 14 gausses/inch/amp = 0. 5 173 x 10‘ 3 tesla s/m/amp 


(IV- 14) 


IV.2 Measured Data 

Coil resistances are measured to check the system characteristics before the 
magnetic field measurement. The results are given in Table IV. 1. The field intensities of 
various coils are measured using a BELL640 Incremental Gaussmeter. For all tests, the 




Table IV. 1 Coil Resistances 


Coil Name 

R(«) 

Helmholtz Coil (Hx - Bias field) 

1.5 

Drag Coil (Hxx -Gradient Field) 

2.1 

Inner Saddle (H^)) 

0.3 

Outer Saddle (H^outj) 

0.5 

Lift Coil (Hzx) 

1.0 

Side Force Coil (Hyx) 

1.0 


Table IV.2 Bias Field Strength at I = 3 amps, y = 0, z = 0 
(Plotted in Figure VI. 8) 


Position x 
(inches) 

-3 

-2 

-1 

0 

1 

2 

3 

Field 

Strength B 
(gausses) 

43 

43 

43 

43 

43 

43 

43 


Table IVJ Gradient Field Strength at I - 3 amps, y = 0, z = 0 
(Plotted in Figure IV. 8) 


Position x 
(inches) 

-3 

-2 

-1 

0 

1 

2 

3 

Field Strength 
B 

(gausses) 

15.6 

10.5 

5.9 

0 

-6 

-10.2 

-15.6 
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current applied to the coil is 3 amps. The direction of the vector B is found by rotating the 
probe until the output reaches a ma ximum 
IV.2.1 Helmholtz Coils 

For the Helmholtz coil system, the bias and gradient fields were tested separately. 
Experimental data are listed in Table IV. 2 and IV.3, and plotted in Figure IV.8. 
Comparison with the computational results given by Equations IV. 4 and IV. 7 shows the 
performance of the Helmholtz coil system is satisfactory. 

B x /I = 14.33 gausses/amp = 1.433xl0' 3 tesla s/amp (IV- 15) 

Bxx /I = 1.67 gausses/inch/amp = 6.57xl0' 3 tesla s/m/ amp (IV- 16) 

IV.2.2 Saddle Coils 

One of the inner saddle coils was tested, and the results were assumed to represent 
the performance of the entire coil system. Experimental data are listed in Table IV. 4 and 
plotted in Figure IV. 9. The magnetic field per amp generated by the half inner saddle coil 
is 

Bt(ia) /I= 3. 1 gausses/amp = 0.3 1 x 10‘ 3 tesla s/amp (IV- 1 7) 

Note that this value is very close to the estimated value given in Equation (IV- 11). 

IV.2.3 Lift and Side Force Coil Assemblies 

Only one of the eight coils which comprise the lift and side force assembly is 
tested. The data are given in Table IV. 5. The performance of a pair of coils and the whole 
lift system (four coils total) are plotted in Figure IV. 10a and b by superimposing the fields 
generated by the other three coils assuming the same characteristics. The slope of Bz in 
the test section is quite uniform, which is shown in detail in Figure IV. 10c. 
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Bxz /I = 0.987 gausses/inch/amp = 3.886xl0' 3 tesla s/m/ amp 


(IV- 18) 


This value corresponds to the estimated value given in Equation (IV- 14). The significant 


discrepancy between this measured value and the estimated value is discussed in Section 


IV.4. 


Table IV.4 Half Inner Saddle Coil Field Strength at I = 3 amps, y = 0, z = 0 

(Plotted in Figure IV. 8) 


x (inches) 

mm 

-2.5 

0 

2.5 

mm 

Field Strength B 
(gausses) 

2.7 

9.3 

9.3 

9.3 

2.25 


Table IV.5 One Lift Coil Field Strength at I = 3 amps, y - 0, z = 0 

(Plotted in Figure IV. 9) 


x (inches) 

-17.5 

-12.5 

-10 

H 

-2.5 

0 

2.5 

H 

12.5 

17.5 

Field Strength 
B (gausses) 

12 

22.2 

24 

19.8 

8.6 

1 

3 

1.14 

0.5 

0.36 

j 

! 


IV.3 Computational Results 

A finite element method of computing the magnetic flux density is developed to 
provide designers with a good reference. This method can also properly represent the iron 
core effect which should explain the discrepancy between the estimated values and 
experimental data. The finite element code used to compute the magnetic fields is 
TOSCA, by Vector Fields Incorporated. This software is a commercial product which has 
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been used by NASA for several years. By inputting the geometries of the coils as well as 
the positions and geometries of the iron cores, the code gives the magnetic flux intensity 
distribution in a three-dimensional space. The three-dimensional stationary 
electromagnetic fields Bx, Bxx, and Bxz are computed. Results for the transverse field B^) 
and Bkoud have not been calculated at this time. Figure IV. 1 1 shows the model of the 
Helmholtz coils as well as the bias field density. Figure IV. 12 gives the gradient field 
created by the anti- symmetric coil windings; note that the variation of the magnetic field in 
the test section (-0.0768 m to 0.0768 m) is approximately linear. The field intensity, per 
amp, of a pair of coils in the lift and side force assembly is computed using a model which 
includes both the coil and the iron core (Figure IV- 13). The model assumes a simplified 
geometry which neglects the effects of fillets in the comers and the copper plating of the 
iron core poles. The computed fields are given in Equations (IV- 19, 20 and 21) below: 


Bx /I = 1 1.78 gausses/amp = 1. 178xl0' 3 tesla s/amp 


Bxx /I = 131 gausses/inch/amp = 5. 15 58x1 O' 3 teslas/m/amp 


Bxz /I = 1.52 gausses/inch/amp = 5.98xl0' 3 teslas/m/amp 


(IV- 19) 
(IV-20) 
(IV-21) 


IV. 4 Discussion 

The three sets of results are compiled in Table IV. 6. Note that the computational 
and test values are all found to be within 44% of the estimated values except for Bxy 
H owever, the computational result, which includes the iron core effect, is very close to the 
test value. Thus one conclusion is that the magnet system still works as it was originally 
designed. Another conclusion is that the iron cores can greatly increase the magnetic field 
strength, in this case about 7.5 times. 




For B x , it is not unusual to find that the estimated value exactly matches the 
computational value, since the simple geometry of the Helmholtz coils allows the flux 
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density to be estimated using the same exact integral which is used to obtain the 
computational value. However, these two sets of data do not take into account the iron 
core effect brought on by the coil assemblies located adjacent to the Helmholtz coils. The 
same effect can be seen in B*x, although the estimated and computational values differ 
slightly due to the simplified geometry used for the estimated value. The relatively large 
discrepancy between the estimated and measured values of B, can be attributed to the 
greatly simplified geometry (using coils of infinit esimal cross-section) assumed for the 
estimation. Note that a simplified geometry, which neglects the influence of the copper 
plating of the iron core poles, is also used to obtain computational values for B*y. This 
assumption could result in the systematic error of the computed value. 


Table 1V.6 Comparison of Three Sets of Results for Magnetic Flux Density 



Estimation 

Measurement 

Computation 

Bx /I (teslas/amp) 

1.178xl0‘ 3 

1.433xl0* 3 

1.178x10° 

Bxx /I (teslas/m/amp) 

5.27x10 3 

6.57xl0‘ 3 

5.1558x10° 

B«in) fl (teslas/amp) 

0.55x10 3 

0.31xl0* 3 


B«out) /I (teslas/amp) 

0.558xl0‘ 3 



Bxz /I (teslas/m/amp) 

0.5173xl0' 3 

3.886xl0* 3 

5.98x10° 
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Two other tests are performed in order to illustrate the effect of the iron cores. 

The first is a measurement of the magnetic flux through the test section created by one left 
coil at I = 3 amps. The result of this test is shown in Figure IV. 14. Note that the 
measurement is identical to the design requirement shown in Figure m.3b. The next test 
involves measurement of the magnetic field strength distribution around the coil with the 
iron core. The results are shown in Figure IV. 15. Again, note the magnitude of the field 
strength on the side toward the tunnel (zero-degree direction) is much higher than on the 
other side due to the iron core effect. 



Figure IV. I Magnetic Field Generated by a Single Loop 



Figure IV.2 Helmholtz Coil Model 









x (inches) 

Figure IV.8 Measurement: Magnetic Field by Helmholtz Coil System at I = 3 A 



x (inches) 


Figure EV.9 Measurement: Magnetic Field by Half Inner Saddle Coil System at I = 3 A 
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Figure IV. 11 Computational Results: Bias Field by Helmholtz Coil System 
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CHAPTER V 

CONTROL SYSTEM FOR THE ASPS 

The Annular Suspension and Pointing System at the Department of Aerospace 
Engineering, Old Dominion University, was originally developed by Sperry Flight Systems 
(now Honeywell Satellite Systems) in the 1970s. Most of the system was loaned to ODU 
in 1992 by NASA Langley Research Center. The original analog control hardware of this 
system was replaced by a new digital controller. The hardware status of the Annular 
Suspension and Pointing System (ASPS) is described in detail by NefF[10] and Neff and 
Britcher [11], Most of the modifications to the system since then have been related to the 
software. A general outline of the system, as well as detailed progress in the controller, is 
given in this chapter. The intention is that the control software be kept somewhat generic 
and used as the basis for the 6-inch MSBS controller in the future. 

V.l System Description 

The ASPS, shown in Figure V.l, consists of a payload mounting plate, a Vernier 
Pointing Assembly (VP A), two coarse gimbal assemblies, a mounting and jettison 
assembly, control electronic rack and balance and testing fixture. The payload mounting 
plate is a machined aluminum structure with a mass of 12.02 kg, a diameter of 965 2 mm 
and thickness of 22.2 mm. It can be mounted directly to the gimbal mounting platform in 
the event that the VPA is not used during the flight. 

The Vernier Pointing Assembly, shown in Figure V.2, is the main component of 
interest here. The VPA includes an annular iron rotor, five Magnetic Bearing Assemblies, 
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a roll motor, seven pairs of proximity sensors, five vernier latches, a roll resolver, rotary 
transformer, and battery packs. The suspended model is an annular iron rotor with a mass 
of 21.6 kg. The system has three axial magnetic bearings, named A, B, and C, spaced 120 
degrees apart to levitate this rotor. The lateral position of the payload is controlled by two 
radial magnetic bearing assemblies, named U and V, which are located 90 degrees apart 
and ±45 degrees from the axial bearings. The roll motor is on the symmetric axis of the 
two radial bearings and opposite to one of the axial bearings. The five vernier latches are 
used to support the iron rotor during launch and recovery maneuvers. 

Each magnetic bearing assembly uses a pair of sensors to detect the rotor position. 
An oscillator provides a 10 MHz, 30 V peak-to-peak signal used by the fourteen proximity 
sensor coils. The proximity sensors used on the ASPS are model KD-2300-12CU, a 
product of Kaman Instrumentation Corporation. The electromagnetic field generated by 
the sensor coils couple with the rotor, and changing gaps cause the impedance of coils to 
vary. The differential signals from each pair of sensors are then sent to a demodulator card 
to realize signal processing. The processing includes amplification, filtering, and 
demodulation and finally gives a linear output voltage. 

It is interesting to compare the position sensors used in the two magnetic 
suspension systems. The proximity sensors used on the ASPS are off-the-shelf products, 
while the EPS used on the 6-inch MSBS is an experimental device developed specifically 
for this application. Also, the ASPS sensors can only be used to detect motion in a single 
degree of freedom, whereas the EPS can resolve motion in five degrees of freedom 
simultaneously. Both types are non-contacting sensors which use the principle of 
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impedance variation, where the coupling between a coil in a sensor and a target is 
dependent upon relative displacement. 

The latches can locate the rotor at the center position in case the MBAs fail. To 
prevent damage to the MBAs, position sensors, and other delicate electronics when not in 
use, the rotor can be locked in the down position. 

The roll resolver and rotary transformer are located at the center of the VPA. 

These two systems are used to measure the angular rotation and position of the payload 
The rotary transformer is used to transfer power between the stationary base of the VPA 
and the payload. These two parts and the roll motor, which complete the roll control, are 
not used at the present time. 

The control electronics rack is a free standing rack that connects to the VPA via 
flexible cables. It contains the power supplies and amplifiers to drive the MBAs and 
sensors. Two Copley Controls Series 661 unregulated D.C. power supplies produce 36 
volts D C. at a rating of 12 A maximum to power the six amplifiers which are needed to 
control the actuators. A voltage of ±15 volts D C. is also needed to drive the oscillator 
card. Copley Controls Corp. 303 series amplifiers are applied in this system. 

The controller, which includes data recording, computation and readout functions, 
is realized by a digital computer, a timer, and two data acquisition boards. Keithley Data 
Acquisition provides two interface boards to realize A/D and D/A conversions. The DAS- 
1402 is a multifunction, high-speed, analog-to-digital interface board for IBM and 
compatible computers. This board offers 8 differential or 16 single-ended analog inputs 
with 12-bit resolution at 100 k samples/sec. The DAS- 1402 can be set for both unipolar 
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(0-10 V) or bipolar (±10 V) signals. DDA-06 is a digital-to-analog I/O expansion board 
providing 6 channels of 12 bit analog output and 24 lines of digital I/O. A timer board is 
also needed to control the sampling rate and handle interrupts. A DT2819 counter/timer 
board includes an Am9513A System Timing Controller, and provides a software- 
selectable 5 MHz or 1 MHz base frequency and a single interrupt which is configurable 
from level 2 to 7. 

V.2 Model Establishment 

The block diagram in Figure V.3 shows the ASPS control scheme. The position 
information of the rotor is acquired by proximity sensors. The signals from the sensors are 
then converted from analog to digital and compared to the commanded position to derive 
the error signals. Input data are processed via a proportional plus derivative controller. 
The output signals are then converted back into analog form to drive the actuators. The 
position of the rotor is sensed continuously and fed back to the computer to maintain 
system stability. 

Each dynamic component in the feedback loop (shown in Figure V 4) has been 
established discretely either through known physical features of the hardware or by testing 
[10, 11], The parameters are listed in Table V. 1. Note that, to simplify the system, the 
transfer functions of the actuators relating the current (i) to the force (F), have been 
linearized near the steady-state point [10, 11], The open-loop transfer function of the 
forward path is: 

349.49 

G(s) = ( K p +K d s ) 00327s 3 + 7.19s 2 (V-l) 



The design process for the PD controller is summarized in Figure V 5 Note that 
due to unit conversions, the actual values of the gains used in the code differ from those 
given in the controller transfer functions. In Controller I, the values of the proportional 
gain Kp and derivative gain Kd were determined by Neff and Britcher [11] using the 

K 

MATLAB Control System Toolbox. Choosing the ratio — — = 0.1 yields 

K p 

. __ 349.49 + 34.949s (V-2) 

(j(S) = K„ - 

p 0.0327s 3 + 7.19s 2 

From the root locus of the system shown in Figure V.6, it can be seen that the closed-loop 
system will remain stable for a range of Kp from 0 to +oo. The values of Kp and Kd were 
chosen to be 4.5 and 0.45, respectively, to give a damping ratio of 0.7. However, the real 
system cannot maintain stability using this design. The actual values used in Neffs code 
[10] were 403 and 5.5 for Kp and Kd, respectively. These values were determined 
experimentally, and the new controller using these values is designated Controller II. 


Table V.l Constants Used in Transfer Functions 


k a 

X 



N 

Io 

s 

go 

1.021 

0.00455 


810 

1.25 

1641.7 

3.41 

A/V 

sec/rtd 


turn/ coil 

A 

mm 2 

mm 

m 

mmsm 

mm m 

Ks(C) 

■Sun 

... Ks(V) 

Ks(W) 


0.00284 

0.00308 

0.00456 

0.00115 

0.00111 

0.00100 

■m 

mm/count 

mm/count 

mm/couflt 

mm/count 

mm/ count 

mm/count 


This work begins with the verification of the transfer functions used by Neff [ 1 0] 
The model parameters are checked and found correct. The next step is a retuning of the 
system using an improved code, which is introduced in the next section. The new code 
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gives values of Kp and Kd of 440 and 4.5, respectively. A new controller is then 
introduced using these values, designated Controller III. Note that these gains are very 
close to those in Controller II. It appears that the ratio Kp / Kd should be 0.01, instead of 
0. 1, to levitate the rotor and maintain system stability. The root locus diagram for 
Controller III is shown in Figure V.7. The system can also maintain stability for a range of 
Kp from 0 to +oo. The responses in position of the theoretical controller (Controller I) and 
the actual working controller (Controller III) to a step input are also given in Figure V.8. 
By comparing these two responses, it can be seen that while Controller III has a shorter 
settling time. Controller I has a significantly smaller overshoot. Although Controller I 
appears in theory to be a better design in some respects, it is based upon simulations for 
linear systems and does not take into account the aerodynamic damping factor, determined 
by the physical features of the rotor (large surface areas). Due to these modelling 
inaccuracies. Controller I does not perform well in reality. 

Finally, an integral term is added to the controller to minimize steady-state errors 
(Controller IV). After running this controller with the integral term turned off and tuning 
the system to satisfy the stability requirement, integral control is activated by a switch in 
the software. A suitable value of Ki is determined experimentally to be 3.7. 

The actual process of tuning the controller involves experimentally determining the 
exact values of these gains which yields satisfactory performance. For the axial MBAs, the 
reference voltage corresponds to a current which will generate enough force to offset 
gravity. K P is chosen to provide stiffness, driving the system back to a desired datum, 
when a disturbance is encountered. The primary role of Kd is to maintain system stability, 
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while also affecting the rise time and overshoot. The value of Kd should be around one 
hundredth Kp; an excessive value could cause system vibration. Once the PD controller is 
working properly, the integral term can be activated to minimize the steady-state error If 
the system is moved away from the balance position due to a steady disturbance, in this 
case a non-conducting weight added to the rotor, the gain Kj performs the function of 
driving the rotor back to the balance position within a small time period. 

V.3 Controller Software 

The controller software is written in Microsoft Quick C. The timer board and the 
associated software driver are supplied by Data Translation, and the data input/output 
interface and its software driver are supplied by Keithley Data Acquisition. Both software 
drivers are C-compatible. 

A flow chart comparison of the original program and the revised program is given 
in Figure V.9. In the original code, the controller was included within the main program 
loop and was executed repeatedly until the exit command was given. Since other functions 
such as keyboard input and data output must also be accommodated within this main loop, 
the sampling rate and controller execution were driven by the time required to perform 
these other functions. The sampling rate was thus chosen to be fairly slow at 100 Hz, and 
data output could be viewed only after turning off the controller since it takes relatively 
long time. Although communication between the D/A, A/D, and timer boards was 
established and the iron rotor was successfully suspended, the system was not sufficiently 
robust to resist a significant artificial disturbance, in this case, a nudge away from the 
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commanded position. In real world application, external disturbances might come from 
vibrations of the base, or disturbance forces on the payload. 

The revised program improves the system’s performance in several ways. First, a 
timer board was used to interrupt execution of the main program at the sampling 
frequency of 200 Hz. When interrupted, program execution jumps to a subroutine which 
acquires analog position data, performs the A/D translation, data calculation, and D/A 
translation. This modification was made so that other functions, such as keyboard 
interrupt and data output which require relatively longer periods of time (e g. archiving 
data to disk), would not affect the sampling rate. When the subroutine is completed, 
execution returns to the main program at the point where it left off. Second, an integral 
term is added to the controller to achieve zero steady-state error. Third, a screen output 
table is designed to show the system status while the program is running. All 1 5 
parameters including reference voltage U„f, the proportional gain Kp , the derivative gain 
Kd and the integral gain Ki for the MBAs can be adjusted in real time, during program 
execution to determine the optimum value of these constants. These real time inputs are 
realized by the keyboard subroutine which executes at the end of the data readout loop. 
Also, the increase in sampling rate from 100 Hz to 200 Hz results in better controller 
resolution and is helpful for the reduction of noise. The last 10 sets of data which includes 
all the controller gain and system performance information are saved in array format and 
written into file DATAOUT.ASC. 

The controller tuning process is greatly simplified since real-time data can be read 
from the screen directly. The axial bearings A, B, and C are tuned before the radial 



bearings U and V. Each axial MBA is tuned individually, with the others deactivated, to 
avoid coupling. When the three axial bearings have been tuned, the rotor is levitated and 
the radial bearings U and V are tuned. 

The test results show that the system’s performance has been improved 
remarkably. The system is now more stable, and due to the increase in sampling rate, is 
less noisy and has a greater resistance to external disturbances. Figures V. 10 and V. 1 1 
show the position errors of the five MBAs as a function of time. From the figures it can be 
seen that the position error of the MBA A is within ±0.2 mm while the original error is 
about ±0.5 mm [10]. The errors of MBAs B and C have also been reduced. Errors in the 
radial directions (MBAs U and V) are slightly less, at about ±0.05 mm. Standard 
deviations (root mean square) for five sets of data (from the 5th second to the 10th 
second) are listed in Table V.2. 

The revised program, as well as a testing program which is used to check if the 
timer board works properly, is given in Appendix C. 


Table V.2 Standard Deviations (STD) of the Position Errors 



A 

B 

C 

U 

V 

STD (mm) 

0.0583 

0.0565 

0.1138 

0.0251 

0.0143 
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Figure V.l The ASPS Component Breakdown 
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Figure V.2 VP A Layout 
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Reference e(t) U(t) i(t) F(t) X(t) 



Figure V.3 System Diagram of the ASPS System 



Figure V.4 Transfer Function of ASPS Control System 
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Figure V.5 Controller Design Comparison 

















Real Axis 


Figure V.7 Root Locus of Controller III (all axes have units of seconds' 1 ) 
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Original 

Program 



Revised 

Program 



Figure V.9 Flow Chart of Control Program (secondary loop in revised program accessed 
via timer board interrupt) 
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CHAPTER VI 
APPLICATIONS OF MSBS 


This chapter discusses the application of MSBS as a modem testing facility for use 
in wind tunnels. Problems in the past which have limited MSBS’ applications are also 
discussed. Further some new technology, which gives the MSBS a promising future, are 
considered. The role of this research project in the overall program is also discussed 
briefly. 

VL1 Test Techniques 

Advanced airplane designs drive the development of new testing techniques. 
Advanced testing techniques include the use of adaptive wall test sections, magnetic 
suspension and balance systems, and non-intrusive laser measurement techniques, all of 
which improve the quality of tests. The ability to simulate the aerodynamics of flight is 
therefore constantly improving. The potential of MSBS is valued by scientists due to 
critical requirements such as data accuracy, elimination of aerodynamic interference, an 
adequate range of angle of attack, dynamic test capabilities, and wind tunnel productivity. 

In the earliest wind tunnels, researchers made corrections for some of the support 
interference problems to ensure data accuracy. However, it is now known that secondary 
errors, which are caused by distortion of the flow, still exist and are unpredictable. Errors 
due to model support systems also increase when large areas of flow separation occur on a 
solid structure. Thus various support stmts were designed with aerodynamic 
considerations in mind in order to reduce gross errors in static aerodynamic coefficients 
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such as dynamic pressure, Mach number, and lift coefficient. However, this approach is 
certainly neither a convenient nor an economical way to solve the problem. 

Testing of fighter models usually requires setting the model at high angle of attack. 
It is difficult to maintain force balance measurement accuracy over a large range of angle 
of attack since a single balance cannot easily accommodate the widely varying loads. High 
angles of attack also generate massive separated wakes which can cause excessive support 
interference, and also decreases in the data accuracy. 

MSBS appears to be able to offer a solution to all of the above problems. 
Additionally, increased experience and improvements in MSBS have led to the capability 
of performing a variety of dynamic tests without the use of complex mechanical devices. 
One of the original motivations to generate a magnetic suspension device was to attain 
more accurate dynamic stability data acquisition [2]. Research on unsteady loads or 
vibration could be done by combining the magnetic suspension system with advanced 
sensors and control systems. Applications of MSBS greatly improve the productivity of a 
wind tunnel by saving time and labor associated with sting and balance set up. Once the 
system is calibrated, several tests can be done consecutively to acquire fast and reliable 
data. Table VI. 6 gives some information about the worldwide activities on MSBS as of 
1996 [21]. 

VL2 Problems Arising from MSBS 

The research issue of MSBS emerged in the 1940’s. In 1957, ONERA in France 
developed the first wind tunnel magnetic suspension device in the world [22], After an 
initial surge of interest, however, research on magnetic suspension technology slowed 
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Table VL1 “Operational** MSBSs, 1996 


Organization 

Apporx. 

Size 1 

Application 

Status 

NASA LaRC 

13 inch 

Low- speed drag, R&D 

Inactive 

ODU 

6 inch 

System R&D 

Recommissioning 

Oxford Univ. 

3 inch 

Hypersonic aerodynamics 

Active 

MAI/TsAGI, Moscow 

15 inch 

System R&D 

Inactive 

NAL, Japan 

4 inch 

System R&D 

Active 

NAL, Japan 

23 inch 

System R&D 

Active 

NCKU, Taiwan 

6 inch 

System R&D 

Active 

NUDT, P R. China 

6 inch 

System R&D 

Active 


1. Square- root of wind tunnel test section cross-sectional area. 


down in the following twenty years. The reasons for its lack of acceptance are various. 
Although magnetic suspension of models in wind tunnels offers several advantages over 
conventional mechanical balance systems, such as no physical support affecting flow 
around the model and capabilities for dynamic testing, a major drawback with magnetic 
balances is the large amount of power required to operate them, especially for large 
models at high dynamic pressures. System stability was poor because the sensor and 
control components were complicated and unreliable at that time. The applications were 
limit ed to research and development, since building a large-scale MSBS was considered 
too expensive. Also, the flow disturbance problems caused by mechanical supports were 
of relatively minor concern at that time. 

Today, innovations in both hardware and software for compensation systems allow 
optimization and simplification of the whole system, making h more acceptable. The 
problems associated with attaining flight conditions in a small scale wind tunnel led to the 
application of MSBS in cryogenic wind tunnels [22], 
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VL3 Cryogenic Technology 

The use of cryogenic technology with magnetic suspension equipment is strongly 
linked to high Reynolds number testing. A cryogenic wind tunnel is used to raise the 
Reynolds number (Re) of an aerodynamic test to more closely match the values in flight. 
With a decrease in temperature, working gas density p increases and viscosity p decreases 
which greatly increases the Reynolds number. Cryogenic operation also lowers the 
velocity for a given Mach number. Taking advantage of MSBS could avoid serious 
interference effects which may mask the effects of variations of Reynolds number. 

The U. S. National Bureau of Standards defines the field of cryogenics to involve 
temperatures below 123K (-150° F). The boiling points of “permanent gases ’’ (helium, 
hydrogen, neon, nitrogen oxygen and air) all lie below 123K. Nitrogen is often used as the 
working gas in cryogenic tunnels [23]. 

VL4 Applications of 6-Inch MSBS 

Following the successful development and demonstration of cryogenic tunnels, 
research activities relating to the 6-inch MSBS sponsored by the NASA Langley Research 
Center were resumed in 1976. During the period from 1976 to 1981, research was 
performed at the Massachusetts Institute of Technology to study the feasibility of using 
superconductors for magnetic balance coils [24]. A superconductor is one of a series of 
metals or compounds which, when cooled to cryogenic temperatures, offers no resistance 
to the flow of electric current. Thus, superconducting magnets have the potential to 
drastically reduce the power requirements of a Magnetic Suspension and Balance System, 
possibly making it feasible to use in larger wind tunnels. NASA LaRC spent the period 
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from January 1977 through August 1980 studying the feasibility of using superconductors 
for a general-purpose magnetic suspension and balance system [24], Several designs and 
specific recommendations were given for the 6-inch MSBS reformation. It soon became 
obvious that the emergence of new technology in the field of superconductivity, large 
magnet construction, computers, control techniques, and innovation in MSBS systems 
made the application of MSBS to large tunnels, not only feasible, but very attractive. 

Based on the previous discussion, further development of the 6-inch MSBS will 
focus on its application to high Reynolds number facilities, either a pressurized tunnel or 
cryogenic tunneL This development would lead to the long-term goal of the integration of 
the 6-inch MSBS into a complete facility. The work which has been done so far involves 
the initial evaluation of the existing MSBS hardware, and development of a digital 
controller for the ASPS system, which can later be modified for use in the MSBS. 
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CHAPTER VH 
CONCLUSIONS 

Two systems, the 6-inch Magnetic Suspension and Balance System and the 
Annular Suspension and Pointing System, have been further developed in support of the 
magnetic suspension research program at ODU. The main focus is on the 6-inch MSBS, 
which is in the first stage of the recommissioning work. 

From the comparison of the estimation, testing, and computation of the magnetic 
field strength distributions of the 6-inch MSBS, it can be concluded that the magnet 
system is still working properly. The inclusion of iron cores in the coil assemblies was able 
to increase the field strength effectively by up to 750 percent. Efforts have been made to 
reinstate the Electromagnetic Position Sensor; however, further documentation and an 
amplifier are needed before the system can be reactivated. The circuit diagrams of the EPS 
have been analyzed and its design has been documented to aid in future work. 

A majority of the improvements to the ASPS have been made in the controller 
software. A timer board has been implemented to enable a higher sampling rate, thus 
improving system performance. The system now operates with less noise and is more 
robust. The timer board also allows more flexibility in performing other tasks such as data 
output and keyboard interruption. The original PD controller has been upgraded to a PID 
controller to minimize the steady state error. Screen output is now clearer and gives a 
direct, real-time readout of system performance. This digital controller has successfully 
replaced the ori ginal analog control system and fulfills the iron rotor suspension and 



97 

balance task. This revised controller is somewhat generic and will be helpful for the new 
design of the 6-inch MSBS control system. 

The future work on the 6-inch MSBS will focus on the further recommissioning of 
the hardware. The EPS system must be calibrated after a proper amplifier is installed. The 
power supply system, cooling system, and control system must also be redesigned to 
accommodate new applications of the 6-inch MSBS, for example in high Reynolds number 
wind tunnels. The final step will be reassembly of the wind tunnel components. 

Future work on the ASPS system should include the reactivation of the sixth (roll) 
degree of freedom, as well as a non-linear controller which can take the non-linear 
characteristics of the actuators into account. Finally, the Balance and Testing Fixture, 
which enables testing in a simulated zero-g environment, can be included and the entire 
system calibrated for actual flight conditions. 
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APPENDIX A 

TRANSFORMATION OF AXES 


Consider two sets of axes shown in Figure II. 1 where a, b, and c are the model’s 
body frame and X, Y, and Z are the wind tunnel frame. The transformation matrix will be 
given by using Euler angles of ‘V” (yaw angle), “0” (pitch angle), and “<J>” (roll angle). 
The axes X, Y, and Z are rotated with respect to the intermediate axes (X', Y', Z' and X", 
Y« 2") following the sequence \j / (about Z), 0 (about Y'), and (about X") to arrive at 
the body axes a, b, and c. 
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APPENDIX B 

EPS CIRCUIT DIAGRAMS 


The complete circuit diagrams of Electromagnetic Position Sensor (EPS) is shown 
in this appendix. Blocks are zoomed in to show details. Selected labels are given to 
indicate certain functions realized by typical elements, such as amplifiers and filters. 
Footnotes are explained at the end of this section. 
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Figure B.l EPS Circuit Diagram 
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Figure B.2 Circuit Diagram, EPS Block 10 
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Figure B.3 Circuit Diagram, EPS Block 1 1 
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Figure B.4 Circuit Diagram, EPS Block 12 
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Figure B.5 Circuit Diagram, EPS Block 13 
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Figure B.6 Circuit Diagram, EPS Block 20 
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Figure B.8 Circuit Diagram, EPS Block 22 
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Figure B.12 Circuit Diagram, EPS Block 32 
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Figure B.16 Circuit Diagram, EPS Block 42 
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Figure B.17 Circuit Diagram, EPS Block 43 
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.19 Circuit Diagram, EPS Block 5 1 
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Figure B.20 Circuit Diagram, EPS Block 52 
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Figure B.22 Circuit Diagram, EPS Block 54 
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Footnotes for Figures B.1-B.22 

1 This is a Bandpass-Low Q-Resistive input filter [25], It is used in serial to eliminate the 


noise from oscillator. 

2. This is a passive high pass filter. It also used for all 7 input signals, (see Block 12, 22, 
32, 42 52 and 53) 

3. This is a typical inverting amplifier. Some other noninverting amplifiers are also 
used [26], 

4. The quadrature oscillator components are produced by active 90 degree phase lag 
circuitry, inverters, and summing networks [17]. Vernier adjustments are also 
available. 

5. The demodulator reference and nulling voltages are generated as sums of in-phase U1 
and quadrature oscillator components U2. 

6. U1 and U2 are added together by a summing circuit to get the demodulator reference 
signal which will go to AD 630 demodulator (See Block 13). 

7. D coarse and Q coarse are adjustable resistors shown on front panel of the electrical 
box They can be tuned to get the nulling signals. 

8. This is also on the front panel of the electrical box. 

9. This is an amplify with adjustable gain show on the panel. 

10. The 2nd adjustable gain can be reached inside the electrical box on each function 
board. 

11. This is the AD 630 synchronous demodulator chip [17], 

12. After signals are demodulated, a low pass filter is added to eliminate high frequency 


noise. 
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13. This adjustable resistor is used to offset the output signal. It can be reached inside the 
electrical box. 

14 Third gain is on the panel too. This is the last and only gain for the D C. output signal 



APPENDIX C 


SOFTWARE 


C.l Program Testl.c 
/* TEST1.C 

This program is developed to test the interrupt ability of 
the Pacer T ime r Board. This program will print out "I am 
doing some less important jobs" 80 times in a certain time 
interval. During this period, every two seconds the time 
board sends out a low pulse signal which drive the excution 
to the interrupt subroutine to print out "Do a more important 
job". 

In order for this program to work the Dual Inline Package 
(DIP) switches on the board surface must be set in a proper 
locations such that the I/O base address is 230 hex (default 
vaule) and the interrupt priority level is at level 4 instead 
of default value of level 5. Also, the output of counter 3 
must be connected to the interrupt input pin for the 
measurement routine to work. 

To compile this program correctly it must be linked with 
PACER. LIB. When compiling a user defined subroutine, stack 
checking must be disabled. 

Written by: Yang Yan 
May 10, 1997 

Annular Suspension and Pointing System 
Old Dominion University 

Department of Aerospace Engineering */ 


#define INTLEVEL 4 

//include <stdio.h> 
#inchide <conio.h> 
//include <math.h> 


//include "pacerdefh" 
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#pragma check_stack(off) 

WORD mode, count; 
double counttime; 
int far k,m,n; 
unsigned long icount; 

void start_measure(int seconds); 

int far impt_proc(void); /* Important job subroutine */ 

main(void) 

{ 

int result; 
double freq; 

PA_ENITIALLZE( 0x230); /* Base address set up */ 

PA_RESET(); 

PASETUPfNTERRUPT(INTLEVEL,impt_proc); 

freq=2000; 

PASETFOUT(freq); /* Set Fout to the desired frequency */ 


start_measure(2); 

print^"all jobs are done"); 

PA_TERMINATE_INTERRUPT(); 

retmn(O); 

} 

void st art mea sure( int seconds) 

{ 

WORD mode,count; 
double freq; 

WORD status; 
unsigned long icount; 


/* inititialze the first counter */ 

PA_BUILD_MODE(MODE_D,SOURCE_F4,GATE_NONE,MODE_OUTPUT_LOW_ 
PULSE, &mode); 

PA_REPORT_FREQUENCY(3 ,&freq); 
count=floor( freq* seconds +0.5); 
counttime=count/freq ; 
if( fabs(coimttime-seconds)>0.2) 



printf("\n measure period to long\n"); 
PA_INmALIZE_COUNTER(3,mode, count, 0); 

PAEN AB LEINTERRUPT ( ) ; 
PA_WRITE_COMMAND(armCounters,Ox4); 


m=0; 

do 

{ 

m=m+l; 

a=0; 

printf("I am doing some less important jobs\n"); 
do 
{ 

n=n+l; 

k=0; 

do 

{ 

k=k+l; 

} 

while (k <= 30000); 

} 

while(n <=11); 

} 

while(m <=79); 


} 

int far impt_proc(void) 

{ 

PA_CLEAR_INTERRUPT(); 

PADIS ABLEINTERRUPT( ); 

PACHECKINTERRUPT(&icount); 

print^"icount=%d, get an interrupt signal\n",icount); 

printf("Do a more important job\n"); 

PA_ENABLE_INTERRUPT(); 

retum(O); 

} 



C.2 Program Testbcd.c 


/♦ CODE: testbcd.c 

This program is developed to use the Pacer Timer Board to perform 
a controled loop within a program A PID controller is designed 
and execute as an f=0.005 Hz interruption. The parameters of the 
controllers can be changed by keyboard interruption. A colorful 
screen output is applied to observe the system performance. The 
final report is written into DATAOUT. ASC which includes the last 
ten sets of parameters of the controllers and system perfomance 
information. 

Controller subroutine: 

This program will read six channels from the DAS- 1401 A/D Data 
acquisition board. The input values will be used to calculated the 
distance, velocity, and integal of the distance of the iron ring. 

Multiply by the gains, the program writes six outputs of controllers 
using the DDA-06 D/A data acquisition board. 

Keyboard subroutine: 

This program will adjust the values of reference voltages to suspend 
the ring, proportional gains, integral gains, and derivitive gains 
in five degrees of freedom Also, a integral controller switch is 
arranged so that it can be turned on after system is started and 
remain stable. 

Note: 

FOUT pin must be connected to SOURCE 1 in order for this program to 
work properly. 

To compile this program correctly it must be linked with several 
libraries which are: DAS1600.LIB, DASRFACE.LIB, PACER.LIB. The file 
myl402.cfg must be in the C:\qc25\bin directory before compiling 
this program This file contains the information about the DAS- 1402 
configuration. 

This program should run from DOS command line. 


Written by Yan Yang 
May 31st, 1997 

Annular Suspension and Pointing SystemD 
Old Do mini on University 
Department of Aerospace Engineering */ 
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#define INTLEVEL 4 

^include <stdio.h> 

#inchxde <conio.h> 

^include <graph.h> 

#inctude <math.h> 

#include <float.h> 

#indude <dos.h> 

#inchide "userprot.h" /* Include file for Keithley DAS- 1402 D/A Board*/ 

#inchide "pacerdef h" /*Inchide file for Pacer Timer Board */ 

/* Note the order of above Include Files */ 
flpra gma check_stack(off) /* This is necessary for interrupt output */ 

^define ESC 27 
#define UP 72 
#define DOWN 80 
#define RIGHT 77 
#define LEFT 75 

#define base 400 /*Base address of DDA-06 D/A Board*/ 


/* Define data type */ 

WORD mode, count; 

double freq,seconds=0.005; /* This is the testing pace */ 
double freq3=50000; /* Frequence to be set at the FOUT pin */ 

/* This is the frequence used by Counter 3 */ 

DDH DAS 1402; /* DevHandle contains base address of DAS 1402 */ 

FRAMEH frameHandle; 

long far clkTicks; 

int high,low,i,counter; 

unsigned result,MODE; 


char ch,chin,NumberOfBoards; 

char buffer[80]; /* Output table array */ 

float far data[10][6][4]; 

char name[ 1 0]= { " ABCUVW" } ; 

int pg; 

char AAA; 

long ADvahieA,RefV ohsA,XdemandA,XoldA,ControlA; 
int highA, lowA; 

float Velocity A,XoffsetAControlVohsAhtVohsA; 



float GainControlDA,GainControlPA,GainControlIA; 


char B,BB; 

long ADvahieB,RefV ohsB,XdemandB,XoldB,ControlB ; 
int highB,lowB; 

float VelocityB,XofisetB,ControlVohsB,InVohsB; 
float GainControlDB,GainControlPB,GainControllB; 

char C,CC; 

long ADvahieC,RefVoksC,XdejnandC,XoldC,ControlC; 
int highC,lowC; 

float VelocityC,XoffsetC,ControlVohsC,InVoltsC; 
float GainControlDC,GainControlPC,GainControlIC; 

char U,UU; 

long ADvahieU,RefVohsU,XdemandU,XoldU,ControlU ; 
int highU.lowU; 

float VelocityU,XoflsetU,ControlVohsU,InVohsU; 
float GainControlDU,GainControlPU,GainControlIU; 

char V,W; 

long ADvahie V,RefVohsV,XdemandV,XoldV,ControlV ; 
int highV,lowV; 

float VelocityV,XofifeetV,ControlVoksV,InVohsV ; 
float GainControlDV,GainControlPV,GainControlIV; 

char W,WW; 
int highW,lowW; 
int ControIW; 

intnMij; 

float dda[2000],ddb[2000],ddc[2000],ddu[2000],ddv[2000] ; 
int o; 

int IntSwitchN); 

float IntegratorA,IntegratorB,IntegratorC,IntegratorU,IntegratorV ; 

char on[ll]="Kion 
char ofl[ll]="Ki oflF:-(”; 


FILE *datao; 
FILE *cond; 


void keyboard(void); 
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int far ipt_proc(void); 

main() 

{ 

/* Note order of initializing boards. Pacer before Keithley */ 

WORD result; 

/* Initialize DT28I9 board */ 

PA_IMTIALIZE(0x230); /* Set Timer Board at base address 230 hex */ 

PA RESET(); /* Set board to start up conditions */ 

PA_SETUP ENTERRUPT(INTLEVEL,ipt_proc); /* User edit routine */ 
PA_SET_FOUT(freq3); /* Set frequence at FOUT pin to freq3 */ 

/* Base f=5M */ 

/* This will be sent to the Counter 3 */ 

PA_BUILD_MODE(MODE_D,SOURCE_SRCl,GATE_NONE,MODE_OUTP 

UT_LOW_P 

ULSE,&mode); 

freq=freq3; 

count=floor(freq*seconds-K). 5); 

PA_INITIALIZE_COUNTERj( 3 , mode, count, count ); 

/* Initialize DAS- 1402 board and reads myl402.cfg and returns Number of Boards */ 
DAS 1 600_DevOpen( "my 1402. cfg",&NumberOfBoards); 

/* Get Device Handle */ 

DAS 1 600_GetDevHandle(0,&DAS 1402); 

/* Initial values of variables */ 

XdemandA= 1 960; 

XoldA=XdemandA; 

RefVohsA=2900; /* Minumum count value to suspend ring at A */ 
GainControlPA=360. 7 ; □ 

GainControlDA=3 . 67 ; 

GainControlIA=3 .00; 

AA=0; /* Input Channel for MBA A */ 

A=5; /* Output Channel for MBA A */ 

XdemandB=1951; 

XoldB=XdemandB ; 

RefVohsB=2872; /* Min umum count value to suspend ring at B */ 
GainControlPB=3 60. 7 ; 

GainControlDB=3 . 47; 

GainControlIB=3.00; 

BB=1; /* Input Channel for MBA B */ 



B=4; 


/* Output Channel for MBA B */ 


XdemandC=2240 ; 

XoldC=XdemandC ; 

RefVohsC=286 1 ; /* Mmumum count value to suspend ring at C */ 

GainControlPC=3 60. 7; 

GainControlDC=2. 97; 

GainControlIC=3 . 00 ; 

CC=2; /* Input Channel for MBA C */ 

C=3; /* Output Channel for MBA C */ 

XdemandU=21 19; 

XoldU=XdemandU ; 

RefVoksU=2507; /* Mmumum coimt value to suspend ring at U */ 
GainControlPU=447. 7 ; 

GainControlDU=4. 47; 

GainControlIU=3 . 00; 

UU=3; /* Input Channel for MBA U */ 

U=2; /* Output Channel for MBA U */ 

XdemandV=1556; 

Xold V=XdemandV ; 

RefVohsV=2373; /* Mmumum count value to suspend ring at V */ 
GainControlPV=447. 7 ; 

GainControlDV=4.47 ; 

GainControlIV=3 . 00; 

W=4; /* Input Channel for MBA V */ 

V= 1 ; /* Output Channel for MBA V */ 


ControIW=2848; 

i=0; 

/* Print out some notes on the screen before it starts */ 
p rintf( "\n\n\n\t\t\tProgr am TESTBCD.C \n"); 
print^ "\t\tAnnular Suspension and Pointing System\n"); 
print^"\t\t\tWritten by: Yan Yang\n"); 
printf("\t\t\t01d Dominion University\n"); 
printf("\t\tDepartment of Aeroqjace EngineeringV); 
printf("\t\t\t MAY, 1997\n\n”); 

print^ "WARNING This program will not run within WINDOWS"); 
prints " only from the DOS command line.\n"); 

printf("\tYou must exit WINDOWS before running this program. \n\n\n"); 
printf("\tConnect the FOUT pin to Source 1 on the Timer board.\n\n"); 



p rintf( "\t\tFrequence set at FOUT pin %.2f\n",freq3); 
printf("\t\tTime to exicute control loop is %1.5f Hz.\n\n", l/freq); 
printf("\t\t\tHit any key to continue\n\n"); 
printff"\t\tHh > to turn on the integral gain.W'); 
printf("\t\tHit < to turn off the integral gain.\n"); 
printf("\t\tThen Hit ESC to end the program.\n”); 
getch(); 


/* Set up screen */ 

_clearscreen( _GCLEARSCREEN ); 
_setbkcolor(6); 

_settextcolor( 1 0); □ 

/* Print the table head */ 

_settextposition( 1,1); 
sprintf(buffer,"NAME H ); 
outtext(buffer); 

_settextposition( 1 , 1 1 ); 
sprintffbuffer/WUTCV)"); 
outtext(buffer); 

_settextposition( 1 ,2 1 ); 
sprintf(buffer,"OUTPUT(V)"); 
outtext(buffer); 

settextposition( 1,31); 
sp rint£(buffer, " dX(mm)" ); 
outtext(buffer); 

_settextposition( 1 ,4 1 ); 

sprinth>uffer,"V(mm/sec)"); 

_outtext(buffer); 

fo r( i=4 j<=\2j.+=2) 

{ 

_settextposition(i, 1 ); 
sprintf(buffer, "%c",name[(i-4)/2] ); 
outtext(buffer); 

} 



_settextposition( 16, 1); 
sp rintf(buffer, "N ame" ) ; 
outtext(buffer); 

_settextposition( 16, 1 1 ); 
sprintf(bufFer, "Ref' ); 
outtext(buffer); 

_settextposhion( 1 6,2 1 ); 
sprint^buflfer,"Kp"); 
outtext(buffer); 

_settextposition( 16,3 1); 
sprint^buflFer/'Kd"); 
outtext(bufFer); 

_settextpoation( 1 6,4 1 ); 

sprintf(buflfer,"Ki"); 

_outtext(buffer); 

for(i= 1 8 ;i<=22 ;i+= 1 ) 

{ 

_settextpositioii(i, 1); 
sprintl^buffer,"%c",name[i- 1 8]); 
outtext(buffer); 

} 


pg=0;D 

o=0; 

PAENABLE_INTERRUPT(); 
PAWRITE_COMMAND(annCounters, 0x4); 


while(ch != ESC) 

{ 

/*print the table*/ 
pg += i; 

if(pg=10) 

{pg-=10;} 


for(m=0 ;m<=5 ;m++) 

{ 



for(n=0;n<=3 31 ++) 

{ 

_settextpositioii(m*2+4,n* 10+11); 
sprintf(bufFer, " %2. 3 f ' , dat a [pg] [m] [n] ) ; 
_outtext(bu£fer); 

} 

} 

_settextposhion( 18, 1 1); 
sprint£(buffer,"%d",RefV ohsA); 
_outtext(buflFer); 

_settextposition( 1 8,2 1 ); 
sprintf(bufFer, "%3 . 3f ',GainControlPA); 
_outtext(buffer); 

_settextposition( 18,31); 

sp rintf(bufFer, " %3 . 3 f ' , GainC ontrolD A) ; 

_outtext(buffer); 

settextposition( 1 8,4 1 ); 
sp rintf(buffer, " %3 . 3f ’,GainControlIA); 
outtext(buffer); 


_settextposition( 19, 11); 

sprintf(buffer,’ , %d ,, ,RefVohsB); 

_outtext(buffer); 

_settextposition( 19,2 1); 
sprintf(bufFer, "%3 . 3f',GainControlPB); 
_outtext(buffer); 

_settextpositioii( 19,3 1); 
sprint^buflFer,"%3 . 3f ' ',GainControlDB); 
_outtext(buffer); 

_settextposition( 1 9,4 1 ); 
sprint^bufifer,"%3.3f',GainControlIB); 
outtext(buffer); 


_settextpositioii(20, 11); 
sprintf(buflfer,"%d",RefVohsC); 
outtext(buffer); 



_settextposition(20,2 1 ); 
sprint£(buffer, "%3 . 3 f GainControlPC); 
outtext(buffer); 

_settextposition(20,3 1); 
sprintf(buffer, "%3 3f \GainControlDC); 
outtext(buffer); 

_settextpositioii(20,4 1 ); 
sprint£(bufiFer,"%3 . 3f ’,GainControlIC); 
_outtext(buffer); 


_settextposkion(2 1,11); 

sprintf(buflfer,"%<l",RefVoltsU); 

_outtext(buffer); 

_settextposition(2 1 ,2 1 ); 

sp rintf(buffer, " %3 . 3f ',GainCoiitrolPU); 

_outtext(buffer); 

settextposition(2 1 ,3 1 ); 
spmtf(buifer,"%3 . 3f ',GainControlDU); 
outtext(buffer); 

_settextposition(2 1 ,4 1 ); 
sprintQbuflfer,"%3 . 3f ’.GainControlIU); 
_outtext(buffer); 


_settextpositioii(22, 11); 
sprintf(buflfer,' , %d",RefVoltsV); 
outtext(buffer); 

_settextpositioii(22,2 1); 

sprint^buflFer, M %3.3f',GamControlPV); 

_outtext(buffer); 

_settextposhion(22,3 1 ); 
sprint£(buflfer,"%3.3f',GainControlDV); 
outtext(buffer); 

_settextposition(22,4 1 ); 
sprint^bu£fer,"%3.3f',GainControirV); 
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outtext(buffer); 

_settextposition(30,21); /* Show integral status */ 
if(IntSwitch=l) 

{ 

sprintf(bufFer,"%s", on); 
outtext(buffer); 

} 

else 

{ 

sprintf(buffer,"%s", off); 
outtext(buffer); 

} 


/* See if get some information from user */ 
keyboard(); 

} 

/* At this point the program returns to the top of the loop 
unless the ESC key has been hit. */ 


fclose(cond); 

PATERMINATEINTERRUPT(); 

PA_RESET(); 

_setvideomode(JDEF AULTMODE); 


for(i=0;i<=5;-H-i) /* Writes a zero volt command to all amplifiers */ 

{ /* to put the amplifies in a no output state. */ 

high=(2047/256); 
low=2047-high*256; 
outp(base+2*i,low); 
outp(base+ l+2*ifiigh); 

} 

/* Final report of all the data */ 

ifi'(datao=fopen("c:\\qc25\\asps\\dataout.asc","w+"))!=NULL) 

{ 

for(i=0^<=9;-H-i) 

{ 

fprintfidatao,"Name\tInV(v)\t\tOutV(v)\t\tX(mm)\t\tV(mm/s)\n"); 

for(j=0j<=4;++j) 

{ 


141 


fprint^datao,"%c\t%2.3f\t\t%2.3f\t\t%2.3f\t\t%2.3f\n",name[j],data[i][j][0],data[i](j][l],d 

ata[i]D3[2],da 

ta[i]D][3]); 

} 

fprintf{ datao, "Wvn 1 '); 

fprintf( datao, "Name\tRefV ( v )\t\tKd\t\tKp\t\tKi\n"); 
fprintf(datao,"A\t%d\t\t",RefVohsA); 
fprintf(datao,"%3.3f\t\t",GainControlPA); 
fprintf( datao, "%3.3f\t\t",GainControlDA); 
fp rintf( datao, "%3.3f\n",GainControlIA); 
fprintf(datao,"B\t%d\t\t",RefVoltsB); 
fprintfl^ datao, ' , %3.3f\t\t , ',GainControlPB); 
fprintf(datao,"%3.3f\t\t",GainControlDB); 

Q)rmt^datao,"%3.3f\n",GainCoiitrolIB); 
fprintf(datao,"C\t%d\t\t",RefVohsC); 
fprintf(datao,"%3 . 3 f\t\t",GainControlPC); 

Q)rint^datao,"%3.3f\t\t",GainControlDC); 

fprintf(datao,"%3.3f\n",GamControlIC); 

fprintf(datao,"U\t%d\t\t",RefVoltsU); 

fprintf(datao,"%3 . 3f\t\t H ,GamControlPU); 

fprintf(datao,"%3 .3f\t\t",GainControlDU); 

fprintf(datao,''%3.3f\n^GainControlIU); 

fpmtf(datao,"V\t%d\t\t",RefVohsV); 

fp rintfi; datao, "%3.3f\t\t",GainControlPV); 

fprintf(datao, ,, %3.3f\t\t",GamControlDV); 

fprintf(datao,"%3.3fji\n",GainControUV); 

} 

fclose( datao); 

} 

else 

{ 

printf("Error:Couldn , t create file.\n”); 

} 

/* NOISE. ASC is the first 10 secs record of the system postion error */ 
cond=fop en( " c : Wqc25\\asp sWnoise. asc", " w+" ) ; 
for( o=0 ;o< 1 999 ;o++) 

Q)rintf(cond,"%3.3f %3.3f %3.3f %3.3f %3.3f\n",dda[o],ddb[o],ddc[o],ddu[o],ddv[o]); 

fclose(cond); 

} 

/* end of the main*/ 


/************** Keyboard subroutine **************************/ 



void keyboard(void) 

< 

if(kbhit() !=0) 

{ 

chin=getch(); 


switch(chin) 

{ 

case ESC:ch=ESC; 
break; 


case'q' : case’Q’: RefV ohsA=RefVohsA+ 1 ; 
break; 

caseW : caseW : RefV ohsA=RefV oltsA- 1 ; 
break; 

case'e 1 : case'E' : GamControlPA=GainControIPA+ 1 ; 
break; 

case , i':case r R':GainControlPA=GaiiiCoiitroIPA- 1 ; 
break; 

caset':caseT:GamControlDA=GainControlDA-H). 1 ; 
break; 

casey :case r Y , :GainControlDA=GainControIDA-0. 1 ; 
break; 

caseY : case' A’: RefV ohsB=RefV oltsB+ 1 ; 
break; 

case's':case'S':RefV ohsB=RefV ohsB- 1 ; 
break; 

case'd': caseD* : GamControlPB=GainControlPB+ 1 ; 
break; 

ca se'f : ca se*F : GainC ontrolPB =GainC ontrolPB- 1 ; 
break; 

case , g , :case'G , :GainControlDB=GamCoiitrolDB+0. 1 ; 
break; 

caseV: case'H" : GainControlDB=GamControlDB-0. 1 ; 
break; 

caseY : case'Z* : RefVohsC=RefV oksC+ 1 ; 
break; 

case'x 1 : case'X 1 : RefV oltsC=RefV ohsC- 1 ; 
break; 

case'c':case'C':GainControlPC=GainControlPC+ 1 ; 
break; 

ca se V : case' V : GainC ontrolPC=GainControlPC- 1 ; 



break; 

easel)' : case'B’: GainControlDC=GainControlDC+0. 1 ; 
break; 

case'll' : case'N 1 : GainControlDC=GainControlDC-0. 1 ; 
break; 

case' 1 ' : case' ! ' : RefVohsU=Re£V ohsU+1; 
break; 

case'2' : case’@': RefV ohsU=RefV oksU- 1 ; 
break; 

case'3' : case'# 1 : GainControlPU=GainControlPU+ 1 ; 
break; 

case'4':case'$':GamControlPU=GamControIPU- 1 ; 
break; 

case'5': case'%' : GainControlDU=GainControlDU+0. 1 ; 
break; 

case'6': case' ^ : GainControlDU=GainControlDU-0, 1 ; 
break; 

case'- : case'+' : RefVolt sV=RefVoItsV+ 1 ; 
break; 

case'-':case'_':Re£V ohsV=RefV oksV- 1 ; 
break; 

case'O': case')': GainControlPV=GamCoiitroIPV+ 1 ; 
break; 

case'9':case'(':GamConttolPV=GainControIPV- 1 ; 
break; 

case'8' : case'*' : GainControIDV=GamCoiitrolDV-H). 1 ; 
break; 

case'T:case'&':GainControIDV=GaiiiCoiitrolDV-0. 1 ; 
break; 

case']' : case' } '; GainControlIA=GamControlIA+0. 1 ; 
break; 

case'f: case' { ' : GainControlIA = GainControlIA-0. 1 ; 
break; 

case , p':case'P:GainControlIB=GainControlIB+0. 1 ; 

break; 

case'o':case'O , :GamControUB=GamCoiitrolIB-0. 1 ; 

break; 

case'i':caseT: GainControUC = GainControlIC+0. 1 ; 
break; 

case*u' ; case'LT : GaiiiCoiitrolIC = CjaiiiCoiitrolIC*0. 1 ; 
break; 

case' ;': case' : GainControlHJ=GainControlIU+0. 1 ; 
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break; 

case'l' : case'L' : GainControUU=GainControlIU-0. 1 ; 
break; 

case'k':case'K':GainControUV=GainControlIV+0. 1 ; 
break; 

case'j , :caseT:GainControlIV=GainControlIV-0. 1 ; 
break; 

case’>' : case' . ' : Int S witch= 1 ; 
break; 

case'^icase'/iIntSwitch^; 

break; 


} 

} 

} 


/**************** interrupt routine ********************************/ 
int far ipt_proc(void) 

{ 

PA_CLEAR_INTERRUPT(); 

PA_DISABLE_INTERRUPT(); 

/* Read the sensor signal: A-D */ 

K_ADRead(DAS 1402,AA,2,&ADvalueA); 

K_ADRead(DAS 1402,BB,2,&ADvahieB); 
K_ADRead(DAS1402,CC,2,&ADvahieC); 

K_ADRead(DAS 1402,UU,2,&ADvalueU); 

K_ADRead(DAS 1402, W,2,&ADvahieV); 

/* Transfer input counts to the voltage value */ 

data[pg][0][0]=((ADvalueA»4)-2047)*20/4096.0; /*InVohsA*/ 
data[pg] [ 1 ] [0]=(( ADvahxeB»4>2047)*20/4096.0; /*InVoksB*/ 

data[pg][2][0]=((ADvalueC»4)-2047)*20/4096.0; /*InVohsC*/ 
data[pg][3][0]=((ADvalueU»4)-2047)*20/4096.0; /*InVoksU*/ 
data[pg][4][0]=((ADvahieV»4)-2047)*20/4096.0; /*InVoltsV*/ 

/* Calculate the displacement, velocity and control values */ 

data[pg][0][2]=((ADvalueA»4)-XdemandA)*.00284; /* XoffsetA in mm */ 
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data[pg][0][3]=(((ADvalueA»4)-XoldA)*.00284)/(couiit/fi'eq); /* VelocityA in 
mm/sec */ 

IntegratorA=IntS\wtch*GainControlIA*data[pg][0][2]*(count/freq)+IntegratorA; 

ControlA=(GainControlPA*data[pg][0][2]+GamControlDA'*data[pg][0][3]+IntegratorA 

+RefVoltsA); 

/*in counts */ 

data[pg][l][2]=((ADvalueB»4)-XdemandB)*. 00308; 

data[pg][l][3]=(((ADvalueB»4)-XoldB)*.00308)/(count/freq); 

IntegratorB=IntSwitch*GainControlIB*data[pg][l][2]*(count/freq)+IntegratorB; 

ControlB=(GainControlPB *data[pg] [ 1 ] [2]+GainControlDB*data[pg] [ 1 ] [3 ]+IntegratorB+ 
RefVohsB); 

data[pg][2][2]=((ADvalueC»4)-XdemandC)*.00419; 

data[pg][2][3]=(((ADvaIueC»4)-XoldC)*.00419)/(count/freq); 

IntegratorC=IntSwitch*GainControlIC*data[pg][2][2]*(count/freq)+IntegratorC; 

ControlC=(GainControlPC*data[pg] [2] [2]+GainControIDC*data[pg] [2] [3]+IntegratorC+ 
RefVohsC); 

data [pg] [3 ] [2]=(( ADvahieU»4)-XdemandU)*. 00115; 

data[pg][3][3]=(((ADvahxeU»4)-XoldU)*.00115y(count/freq); 

IntegratorU=IntSwftch*GainControlIU*data[pg][3][2]*(count/freq)+IntegratorU; 

ControlU=(GainControIPU*data[pg][3][2]+GamControlDU*data[pg][3][3]+IntegratorU 

+Re£VoltsU); 

data[pg] [4] [2]=((ADvalueV»4)-XdemandV)*.00 1 1 1; 

dat a[pg] [4] [3]=((( ADvahieV»4)-Xold V)*. 0011 l)/(count/freq); 

IntegratorV=IntSwitch*GainControlIV*data(pg][4][2]*(countyfreq)+IntegratorV; 

ControlV=(GamControIPV*data[pg][4][2]+GamControlDV*<kta[pg][4][3]+IntegratorV 

+RefVohsV); 


if(o<=1999) 

{ 

dda[o]=data[pg] [0] [2] ; 
ddb[o]=data[pg] [ 1] [2] ; 
ddc[o]=data[pg] [2] [2] ; 
ddu[o]=data[pg] [3 ] [2] ; 
ddv[o]=data[pg][4][2]; 
o=o+ 1 ; 
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} 

else 

{ 

dda[ 1999]=data[pg][0][2]; 
ddb[1999]=data[pg][l][2] 
ddc[ 1 999]=data[pg] [2] [2] ; 
ddu[ 1999]=data[pg][3][2] 
ddv[1999]=data[pg][4][2] 
> 


/"■Convert the output counts to the volts */ 

data[pg][0][l]=(ControlA-2047)*5/4096.0; 
data[pg][ 1][ l]=(ControlB-2047)*5/4096.0; 
data[pg][2][l]=(ControlC-2047)*5/4096.0; 
data[pg][3][l]=(ControlU-2047)*5/4096.0; 
dataQ>g][4][l]=(ControlV-2047)*5/4096.0; 

/* Write the output to the amplify: D-A */ 
highA=(ControlA/2 5 6); 
lowA=ControlA-high*256; 

outp(base+2*A,lowA); /* Writes control output to Amplifier A */ 
outp(baseH+2*A,highA); /* Using DDA-06 D/A Data acquisition Board*/ 
XoldA=ADvahieA»4; /* saves old distance for velocity */ 


highB=(ControlB/256); /* same as above */ 

lowB=ControlB-highB*256; 

outp(base+2*B,lowB); 

outp(baseH+2*B,highB); 

XoldB=ADvahieB»4; 

highC=(ControlC/256); 
lowC=ControlC-highC*256; 
outp(base+2*C,lowC); 
outp(baseH+2*C,highC); 
XoldC=ADvahieC»4 ; 

highU=(ControlU/256); 
lo wU=C ontrolU-highU* 256; 
outp(base+2*U,lowU); 
outp(base+ 1 +2*U,highU); 
XoldU=ADvahieU»4; 


highV=(ControlV/256); 



lowV=ControlV-highV*256; 
outp(base+2*V,lowV); 
outp(base+ 1 +2* VJiigh V); 
XoldV=ADvahieV»4; 

highW=(ControIW/256); 
lowW=ControIW-highW*256; 
outp(base+2*W,lowW); 
outp(base+ 1 +2 * W JiighW); □ 


PAENABLEINTERRUPT(); 

retum(O); 

} 
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